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SECTION  1 
INTRODUCTION 


The  estimating  scheme  described  in  this  report  is  the  direct 
result  of  prior  work  by  Hutchinson,  and  co-workers  '*'-2  and  by 
Shih  and  co-workers  ^  This  estimating  scheme  is  used  to 

determine  a  fracture  mechanics  parameter  (the  J  Integral)  which 
describes  the  intensity  of  the  Hutchinson  ,  Rice  and  Rosengren 
(HRR)  stress-strain  field  at  the  crack  tip  in  an  elastic-plastic 
material.  The  estimating  scheme  can  also  be  used  to  determine 
load-load  line  displacement  and  crack  mouth  opening  displacement. 
The  J-Integral  parameter  determined  by  the  estimating  scheme  can 
be  used  to  obtain  the  crack  tip  opening  displacement  when  the  HRR 
field  singularity  describes  the  crack  tip  singular  behavior. 

The  two-fold  purpose  of  the  report  is  (a)  to  review  the 
theoretical  basis  for  the  estimation  scheme  and  (b)  to  describe 
a  computer  program  which  utilizes  the  estimating  scheme  to  yield 
values  for  the  J-Integral  (J) ,  the  load-line  displacement  (A) , 
and  the  crack  mouth  opening  displacement  (CMOD  or  <$)  for  five 
different  structural  crack  geometries. 
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SECTION  2 


THEORETICAL  BACKGROUND 

This  section  of  the  report  has  been  prepared  to  provide 
background  information  relative  to  the  elastic-plastic  fracture 
mechanics  parameters  which  are  derived  by  the  computer  program 
described  in  Section  3. 

2.1  ELASTIC-PLASTIC  FIELD  PARAMETERS 

The  solutions  for  the  in-plane  tensile  opening  mode  (Mode  1) 

type  of  deformation  problems  are  the  primary  interest  in  the 

present  work  as  well  as  fracture  mechanics  in  general.  However, 

the  mathematical  difficulties  have  prevented  detailed  treatment  of 

elastic-plastic  problems.  Except  in  some  out-of-plane ,  tearing 

9  10 

mode  (Mode  3)  type  of  problems  '  ,  rigorous  mathematical  solu¬ 

tions  of  elastic-plastic  problems  are  not  available  in  general; 
the  available  limited  cases  of  elastic-plastic  crack  tip  stress 
analyses  will  be  reviewed  in  this  subsection. 

Figure  1  shows  the  crack  tip  and  area  ahead  of  the  crack  tip^. 
The  region  ahead  of  the  crack-tip  is  divided  into  three  distinct 
zones:  (1)  elastic,  (2)  elastic-plastic,  and  (3)  intensity  non¬ 

linear  (large  strains  and  rotations,  and  ductile  cavities)  zone. 

The  elastic  zone  (1)  controls  the  behavior  when  the  plastic  zone 
size  is  small  compared  to  the  elastic  zone  and  the  geometry.  In 
this  case,  referred  to  as  small  scale  yielding,  linear  elastic 
fracture  mechanics  (LEFM)  is  applicable.  If  the  plastic  zone 
size  is  large,  compared  to  the  case  of  small  scale  yielding,  LEFM 
is  not  applicable. 

The  intense  elastic-plastic  stress-strain  field  contained 

within  zone  2  of  Figure  1  is  further  expanded  in  Figure  2.  When 

the  intensely  deformed  process  zone  is  small  compared  to  the  size 

of  the  elastic-plastic  zone  under  consideration,  the  deformation 

theory  of  plasticity  for  a  power  hardening  material  can  be  used  to 

obtain  stress-strain  solutions  ahead  of  the  crack  tip  outside  the 

7 

intensely  deformed  process  zone  as  suggested  by  Hutchinson  and 

g 

Rice  and  Rosengren  .  These  authors  expressed  power  hardening 
using  a  stress  (a) -strain  (e)  relationship  given  by: 
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ELASTIC 


/ 


\ 


Zone  1  =  An  Elastic  Field  Surrounding  the  Crack  Tip 

Zone  2  =  An  Elastic-Plastic  Field  Surrounding  the  Crack  Tip 

Zone  3  =  An  Intense  Zone  of  Deformation 


Figure  1.  Crack-Tip  Stress  and  Strain  Fields  Surrounding  the 
Crack  (Reference  11) . 


Figure  2.  Expanded  View  of  the  Elastic-Plastic  Stress-Strain 
Field  (Reference  11)  . 
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where  o  and  e  are  reference  stresses  and  strains,  respectively, 

o  o 

and  N  is  the  strain  hardening  exponent. 

When  Equation  1  is  used  to  model  the  behavior  of  the  material 
in  the  plastic  range,  the  stress,  strain,  and  displacement  functions 
for  the  crack  tip  region  are  given  by  Equation  2,  regardless  of  the 
amount  of  plastic  deformation: 


N 


a.  .  =  K 
iD  a 

a .  . 
ID 

(0,N) *  r 

N+l 

(2a) 

1 

e .  .  =  K 
iD  e 

~£ij 

(0 ,N) • r 

N+l 

(2b) 

1 

u.  =  K 
i  e 

ui  (6, 

XTN  N+l 

N)  •  r 

(2c) 

where  and  K£  are  the  magnitudes  of  singularities  of  appropriate 


quantities  with  K£  being  a  function  of 


The  functions  ck  ^ , 
and  vL  depend  on  angle  and  exponent  N  in  Equation  1.  Equations 
2a,  2b  and  2c  have  been  referred  to  as  the  set  of  "HRR"  field 
equations  after  the  initial  investigators  (References  7  and  8). 


In  the  derivation  of  the  HRR  field  equations,  the  J deforma¬ 
tion  theory  of  plasticity  was  used  to  describe  the  material  behavior. 
If  the  loading  is  proportional,  the  field  solutions  obtained  using 
both  the  ^-deformation  theory  and  the  more  realistic  ^-incremented 
theory  of  plasticity  are  the  same.  However,  when  unloading  occurs 
during  deformation,  the  loading  path  may  be  different  from  the 
assumed  proportional  loading  and  the  validity  of  the  above  solutions 
is  not  guaranteed. 

The  magnitude  of  the  singularity  in  Equation  2  can  be  written 
in  terms  of  J  as  when  the  process  zone  is  small  (Reference  11) : 
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In  equation  4,  r  is  any  contour  that  surrounds  the  crack  tip, 
ds  and  nj  are  the  element  length  and  outward  normal  log  r,  Uj^i 
are  the  displacement  gradients  and  ck  ^  and  e^j  are  the  stress  and 
strain  tensors.  The  J  has  been  determined  to  be  independent  of 
the  location  of  contour  curve  T. 

Hence,  for  a  given  cracked  material,  with  the  assumption  of 
power  law  hardening  behavior,  deformation  theory  of  plasticity, 
and  proportional  loading,  there  exists  a  unique  elastic-plastic 
stress  and  strain  field  which  is  characterized  by  its  intensity  J. 

In  a  manner  similar  to  the  LEFM  approach  where  the  stress  intensity 
factor  (K)  measures  the  intensity  of  stress  and  strain  within  the 
elastic  crack  tip  field,  the  parameter  J  defines  the  intensity  of 
the  elastic-plastic  crack  tip  field,  the  parameter  J  defines  the 
intensity  of  the  elastic-plastic  stress  and  strain  in  the  crack  tip 
field  and  thus  provides  a  basis  for  a  nonlinear  fracture  mechanics 
approach.  The  use  of  J  to  define  the  level  of  elastic-plastic 
stresses  and  strains  around  the  crack  tip  requires  that  the  intense¬ 
ly  deformed  process  zone  is  small. 

For  ideally  plastic  materials,  the  relationship  between  the 
J- integral  and  the  crack  tip  opening  displacement  (6t)  defined  by 
the  opening  distance  between  the  intercepts  of  two  45°  lines  drawn 
back  from  the  crack  tip,  with  the  deformed  profile  of  the  stationary 
crack  as  shown  in  Figure  3  is  given  by  References  12  and  13: 
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Figure  3.  Definition  of  Crack  Tip  Opening  Displacement  6. 
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(5) 


where  d  is  a  constant.  For  strain  hardening  materials,  Shih  and 
n 

co-workers  (References  4  and  14)  have  shown  that  Equation  5  can  be 

used  to  relate  J  and  when  the  constant  dn  is  replaced  with  a 

function  which  is  strongly  dependent  on  the  strain  hardening  exponent 

and  mildly  dependent  on  the  ratio  o0/E  (E  =  Young's  modulus).  The 

functional  relation  for  a  strain-hardening  material  is  described  in 

Figure  4;  Figure  4a  applies  for  plane  strain  conditions  while  Figure 

4b  applies  for  plane  stress  conditions.  Figure  4  has  been  derived 

★ 

assuming  that  Equation  1  describes  the  material.  On  the  basis  of 
Equations  3  and  5,  6t,  the  crack  tip  opening  displacement,  is  also 
a  parameter  which  characterizes  the  intensity  of  the  stress-strain 
HRR  field.  The  above  HRR  field  equations  are  applicable  only  for 
the  case  of  stationary  cracks. 

2.2  PARAMETER  DETERMINATION 

For  elastic-plastic  materials,  the  parameters:  the  J-integral 
(j) ,  the  crack  mouth  opening  displacement  (CM0D=6) ,  and  the  load- 
line  displacement  due  to  the  presence  of  a  crack  (A  )  can  be  approxi¬ 
mated  by  the  contributions  due  to  their  linearly  elastic  and  plastic 
parts  (References  3,  4,  and  16): 


J 

6 


Je  +  Jp 

6e  +  6P 


Ae  +  AP 
c  c 


(6) 


where  superscripts  e  and  p  denotes  elastic  and  plastic, 


respectively. 


* 

If  the  material  is  described  by  the  relationship: 

e  ,  a  .  n  , 

- —  =  a  (-= — )  and  a  =1 

eo  o 

then  dR  is  as  presented  in  Figure  4.  However,  if  a  t4  1  then  d  is 
equal  to  the  product  of  the  value  of  dn  given  in  Figure  4  and  the 
quantity  of  a 1//n  (Reference  15). 
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1.0 


Figure  4. 


(b) 


Variation  of  d  with  N  and  aQ/E  (a)  Plane  Strain 
(b)  Plane  Stress  (Assuming  Equation  1  Applies) . 
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Reasonable  estimates  of  the  total  load-line  displacement  (A)  for  the 
structural  geometry  can  be  obtained  by  summing  the  contributions  due 
to  the  presence  of  the  crack,  i.e.  A^  from  Equation  6,  and  to  that 
due  to  the  structural  geometry  without  a  crack  (An) ,  (Reference  3) : 


A=Ac  +  An  (7) 

The  error  in  using  Equation  7  is  small  when  the  distance  between 
the  load  points  is  much  smaller  than  other  structural  dimensions. 

The  purpose  of  this  subsection  is  to  outline  the  analysis  that 
must  be  accomplished  in  conjunction  with  Equation  6  to  derive 
estimates  of  the  elastic-plastic  parameters.  The  following  three 
paragraphs  present:  the  elastic  formulation,  the  plastic  (strain¬ 
hardening)  formulation  and  the  elastic-plastic  transition  formula¬ 
tions,  respectively. 

2.2.1  Linearly  Elastic  Contribution 

For  linearly  elastic  materials,  the  elastic  crack 
parameters  appearing  in  Equation  6  can  be  expressed  in  the  form: 


Je 


a  e  a 
o  o 


,  o°°  .  2  C.e ,  ... 

=  (— — )  J  (a/b) 


e  a 
o 


=  6e  (a/b) 


(8) 


e  a 
o 


(a/b) 


where  a°°  is  the  remotely  applied  stress,  and  and  £q  are  reference 
stresses  and  strains  related  by  the  expression  a  =  Ee  .  Functions 
J  ,  6,  and  A^  are  functions  only  of  the  ratio  of  crack  length  to 
width  (a/b) .  These  functions  can  be  found  in  Reference  17  for 
different  finite  width  crack  geometries. 
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2.2.2  Plastic  (Strain  Hardening)  Contribution 


To  derive  the  plastic  crack  parameters  given  in 
Equation  6,  certain  assumptions  are  required.  First,  the  material 
is  assumed  to  behave  according  to  a  power  hardening  constitutive 
(cr-e)  equation  of  the  form 


(9) 


where  a  is  a  dimensionless  constant,  a  and  e  are  reference 

o  o 

stresses  and  strains,  and  n  is  the  hardening  exponent.  Note  that  the 
exponent  in  Equation  9  is  in  the  inverse  of  the  exponent  in  Equation 
1.  For  n=l,  the  material  behaves  as  a  linearly  elastic  material; 
as  n  approaches  infinity,  the  material  behaves  more  and  more  like 
a  perfectly  plastic  material.  Generalizing  Equation  9  to  multi- 
axial  states  via  the  J2  deformation  theory  results  in: 


(10) 


where  e. .,  S. .,  and  a  are  the  strain  deviator,  the  stress 
iD  ID  e  _ _ 

deviator,  and  the  effective  stress  (  =  /3/2STTTsTT) ,  respectively. 
Ilyushin  (Reference  8)  first  noted  that  the  boundary  value  problems 
which  (1)  have  an  externally  applied,  monotonically  increasing  load 
or  displacement  and  (2)  are  based  on  Equation  10,  have  some  special 
properties.  He  showed  that  the  stress  at  each  point  in  the  body 
varies  linearly  with  a  single  load  parameter  when  tractions  are 
prescribed  on  all  boundaries  and  the  directions  of  these  tractions 
remain  fixed  while  their  magnitudes  are  everywhere  linearly  propor¬ 
tional  to  the  load  parameter.  Since  the  stress  components  at  each 
point  are  proportional,  the  solution  based  on  deformation  theory 
also  applies  for  incremental  plasticity  theory  when  the  load 
parameter  is  monotonically  increasing. 
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The  functional  dependence  of  the  field  parameters 
(stress,  strain,  and  displacement)  on  the  applied  load  (or  dis¬ 
placement)  also  means  that  crack  tip  parameters  can  be  uniquely 
related  to  the  remotely  applied  load  (a°°)  via  the  following  expres¬ 
sions  (References  3  and  4) : 


a  e  a 
o  o 


■  <S=> 

o 


n+1 


t  ia  \ 

J  (j^ » n) 


6 


P 


£oa 


a  oo 

F 


n 


(11) 


eoa 


-  O 
0 


n 


Ac  (F'n) 


where  Jp,  6P,  and  AP  are  functions  only  of  ^  and  n.  The  reader 
will  note  that  the  functional  forms  given  by  Equation  11  are 
similar  to  Equation  8  and,  in  fact,  they  reduce  to  Equation  8  when 
n=l. 


An  alternate  form  of  Equation  11  that  has  been  used 
previously  (References  4  and  5)  and  which  is  used  in  the  computer 
program  discussed  in  Section  3  is: 


Jp  =  a  aQ  eQ  a  fx  (a/b)  •  h1  (^,n)  • 

6P  =  a  eQ  a  f 2 (a/b)  h2(|,n)  •  (|-)n 

o 

„  ,  p  n 

Ap  =  a  eQ  a  f3(a/b)  h3(^,n)  •  (|-) 


(12) 


where  P  is  the  applied  load  (per  unit  thickness)  ,  PQ  is  the  l-imit  load 
(per  unit  thickness),  f^,f2  an<^  ^3  are  functi°ns  only  of  geometry  and 
crack  length  while  h^,  h2  and  h3  depend  on  geometry,  crack  length, 
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and  strain  hardening  exponent.  Shih  and  co-workers  (References  4 
and  5)  have  tabulated  the  functions  h^ ,  h h^  for  a  number  of 
geometries  for  the  conditions  of  plane  stress  and  plane  strain. 

From  the  reference  tabulated  data,  these  functions  can  be  obtained 
by  interpolation  for  any  value  within  the  ^  and  n  limits  given; 
thus,  the  plastic  (strain-hardening)  component  of  Equation  6  can 
be  computed  for  any  given  applied  load  P  (or  a°°)  from  Equation  12 
(or  11) . 

2.2.3  Elastic  to  Plastic  Transition 

In  the  two  preceding  paragraphs,  we  reviewed  the 
procedures  for  estimating  the  parameters  of  Equation  6  when  the 
material  was  either  linearly  elastic  or  plastically  strain  harden¬ 
ing.  Under  the  assumptions  for  small  scale  yielding,  Bucci,  et  al. 
(Reference  16)  observed  that  the  transition  region  between  fully 
linear  elastic  and  fully  large-scale  plastic  deformation  could  be 
more  accurately  modeled  if  the  physical  crack  length  term,i.e.  "a" 
in  the  elastic  components  of  Equation  6  was  replaced  with  an 
effective  crack  length  term  (ag) .  Bucci,  et  al.  suggested  an 
effective  crack  length  based  on  adding  the  Irwin  plastic  zone  size 
(ry)  to  the  physical  crack  length: 


a 


e 


a  +  r 


Y 


(13) 


where 


(14) 


with  6=2  for  plane  stress  and  6=6  for  plane  strain.  The  stress 
intensity  factor  is  represented  by  K. 

Subsequently,  Shih,  et  al .  (Reference  4)  defined  the 
effective  crack  length  with  a  modification  to  account  for  strain 
hardening  and  for  the  relationship  of  the  load  to  limit  load.  The 
Shih,  et  al.  formulation  for  the  effective  crack  length  is: 


a_  =  a  +  <)>  r 
e  y 


(15) 


12 


where 


2 


(16) 


and 


r 


y 


JL 

Bit 


n-1 

n+1 


<t>  = 


1 _ 

1  +  (P/PQ)2 


(17) 


The  factor  provides  a  correction  to  r^  in  the  case 
of  small  scale  yielding,  and  also  limits  the  contribution  from  the 
plastically  adjusted  crack  length  for  large-scale  yielding  so  that  the 


j  5  A 

values  of  —  ,  —  ,  and  g  approach  unity. 

JP  6P  AP 

c 

The  version  of  Equation  6  embedded  in  the  computer 
program  discussed  in  Section  3  is  given  by: 


J  = 


6 

A 


c 


(ae>  +  jP 

(ae)  +  SP 


(a ,  n) 


(a  ,n) 


(a,n) 


(18) 


where  ag  is  the  effective  crack  length  defined  in  Equation  15.  The 
elastic  contribution  given  by  Equation  8  is  thus  modified  by  the 
replacement  of  the  physical  crack  length  by  the  effective  crack 
length.  The  plastic  contribution  is  given  by  Equation  12. 
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2.3  FORMULATIONS  FOR  DIFFERENT  GEOMETRIES 


This  subsection  has  been  designed  to  present  all  the  rele¬ 
vant  information  and  equations  that  would  allow  one  to  use  the 
estimating  procedure  outlined  in  subsection  2.2  for  five  different 
structural  geometries.  The  equations  and  data  presented  herein 
have  been  previously  presented  by  others  (References  2  through  5 
and  17)  and  are  only  repeated  here  for  completeness.  The  infor¬ 
mation  for  each  structural  geometry  is  summarized  using  a  figure 
and  two  supporting  tables.  The  five  structural  geometries  con¬ 
sidered  are  listed  in  Table  1  along  with  the  individual  figure 
and  table  numbers  which  present  the  supporting  information. 

Each  of  the  five  supporting  figures  all  are  arranged  in  the 
same  way  so  that  each  geometry's  relevant  equations  and  data 
appear  in  the  same  location  in  the  figures.  The  arrangement  of 
the  figures  is  such  that  a  schematic  of  the  geometry  appears  in 
the  upper  left  hand  side  and  the  individual  versions  of  Equation 
18  (the  elastic-plastic  parameters)  can  be  found  in  the  upper 
right  hand  side.  The  remainder  of  the  figure  is  devoted  to 
presenting  equations  that  further  define  the  elastic  and  plastic 
components  of  Equation  18  and  to  defining  the  tables  that  contain 
discrete  value  information  on  the  plastic  functions  h^,  ,  and 

h3  * 

In  order  to  ensure  that  the  information  in  Figures  5  through 

9  is  properly  interpolated,  a  few  additional  notes  are  presented 

here.  First,  E'  is  Young's  modulus  (E)  for  plane  stress  and  is 

o 

equal  to  E/(l  -  v  )  for  plane  strain  where  v  is  the  Poisson's 

T 

ratio  of  the  material.  The  loads,  P  and  P  ,  are  the  applied  load 
and  the  theoretical  limit  load  for  a  perfectly  plastic  material 
(  n  =  °°)  respectively;  both  are  expressed  per  unit  thickness. 
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GEOMETRY 


ELASTIC-PLASTIC  PARAMETERS 


j  =  — y  F*  (z)  |!  +  oo  e  a(L~a)  h,(a/b,n)  (P/P*)n+1 

4b2  E  o  o  b  1  o 

5  =  (z)  +  oteoa  h2(a/b,n)  (P/P^)11 


a  - 
Ac  b 


T,n 


v2  (?)  +aeQ  a  h3<a/b,n)  (P/PQ) 


/TP  ,n 
ouy1E5„ 


An  =  Is  +  73  ae"L(^  1 
A  =  An  +  Ac 


PANEL  (CCP) 

SUPPORTING  ELASTIC  FUNCTIONS:  F,  V3  AND  Vj  with  z  =  ae/b 


F(Z)  =  * 

y  sec  ( 

7TZ  . 

2 

V3(z)  = 

-0.071 

-  0.535z 

+0.169Z2  + 

0.02z3  - 

1 . 071 (1/z ) 

ln(l-z) 

V2(z)  = 

-1.071 

-1.071 

+  0.250z 

(1/z)  In 

-  0.357z2+ 
(1-Z) 

0. 121z3  - 

4 

0.047z  +0. 

5 

008z 

T 

THEORETICAL  PLASTIC  LIMIT  LOAD:  PQ 

Plane  Stress  Condition 
T  ^ 

Po  =/3  ao(b'a) 

Plane  Strain  Condition 

PT  =  2a  (b-a) 

o  o 

SUPPORTING  PLASTIC  FUNCTIONS:  ,  h2,  and  h3 

Plane  stress  condition  tabularized  in  Table  2a 
Plane  strain  condition  tabularized  in  Table  2b 


Figure  5. 


Elastic-Plastic  Parameters  for  the  Center-Cracked 
Panel  (References  4,  5  and  17). 
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and  h3  FOR  THE  PLANE  STRESS  CCP  IN  TENSION  (REFERENCE  5) . 
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AND  h..  FOR  THE  PLANE  STRAIN  CCP  IN  TENSION  (REFERENCE  5). 
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GEOMETRY 


ELASTIC-PLASTIC  PARAMETERS 


J  -  £T  •  p2(z)gf  +  0LOo£o  (b-a)h1(a/b,n)  (P/P^) 


Tx  n+1 


T,n 


6  =  V1(Z)E'  +  aeoa  h2(a/b'n) (p/Po) 


Ac  =  v2(z)l'  +  aeoa  h3(a/b,n) (P/P^)n 


=  0 


COMPACT  TENSION 
SPECIMEN  (CT) 


An  + 

n  c 


SUPPORTING  ELASTIC  FUNCTIONS:  F,  V.^  and  V2  with  z  ■ 


F(z) 


(2+z) 


/z (1-z) 


Y75  (0. 886+4. 64z-13.32z2+14.72z3-5.6z4) 


Vj_(z)  *  (5. 435+43. 315z-83.166z2+57.694z3)/(l-z2) 
and 

V2(z)  =  (0. 995+27. 977z-27.209z2+11.062z3)/(1-z2) 

THEORETICAL  PLASTIC  LIMIT  LOAD:  PT 

o 

Plane  Stress  Condition 
PT  =  1.072n  (b-a) a 

O  0 

Plane  Strain  Condition 

_T  2,52  /w  x 
P^  -  - n  (b-a)  a 

°  /j  0 


with  n  = 

SUPPORTING  PLASTIC  FUNCTIONS:  hj.,  hj,  h3 


^)2  +  2  <&>  +  2  -  +  1) 


Plane  stress  condition  summarized  in  Table  3a 
Plane  strain  condition  summarized  in  Table  3b 


Figure  6.  Elastic-Plastic  Parameters  for  the  Compact 
Tension  Specimen  (References  4,  5  and  17). 
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TABLE  2b 

AND  h  FOR  THE  PLANE  STRAIN  CCP  IN  TENSION  (REFERENCE  5) . 
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FOR  THE  PLANE  STRAIN  COMPACT  TENSION  SPECIMEN  (REFERENCE  5). 
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Figure 


GEOMETRY 


ELASTIC-PLASTIC  PARAMETERS 


J  =  ^  F2(z)|!  +  ao0e0(b-a)h1(a/b,n)  <P/P*)n+1 
5  =  ^  vi^>|'  +  oteo(b-a)h2(a/b,n)  (P/P^)n 

Ac  =  ^  v2  (z) +  c‘e0(b'a)h3<a/b'n>  (p/Po>n 

.  PL  ,  r?  T  ,  /3P  x  n 
n  Eb  o  4bo 


DOUBLE  EDGE 
CRACKED  PANEL  (DCP) 


A  A  +  A 
n  c 


SUPPORTING  ELASTIC  FUNCTIONS:  F,  V^,  and  V2  with  z  =  y 

F  ( z )  =  (1.122  -  0.561z  -  0 . 205z 2  +  0.471z  3  -  0.190 4  )//Irz 


Vi(z)  =  (—■)  [  0 . 459  (Sin  y)  -  0.065  (Sin  y)3  -  0.007  (Sin  y)5 


+  Cosh"1  (Sec  y)  ] 

V2(z)  =  (JL) [0.0629  -  0.0610  (Cos  y)4  -  0.0019  (Cos  y)8 
+  in  ( Sec  y-)  ] 

THEORETICAL  PLASTIC  LIMIT  LOAD:  P^ 


Plane  Stress  Condition 

Po  =  (b'a) 

°  /I  0 


Plane  Strain  Condition 

PT  =  5.94  o  (b-a) 
o  o 

SUPPORTING  PLASTIC  FUNCTIONS:  h^  h2,  and  h3 

Plane  stress  condition  tabulated  in  Table  4a 
Plane  strain  condition  tabulated  in  Table  4b 


Elastic-Plastic  Parameters  for  Double-Edge 
Cracked  Panel  (References  4,  5  and  17). 
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AND  h..  FOR  THE  PLANE  STRESS  DECP  IN  TENSION  (REFERENCE  5)  . 
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AND  h-,  FOR  THE  PLANE  STRAIN  DECP  IN  TENSION  (REFERENCE  5)  . 
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GEOMETRY 


ELASTIC-PLASTIC  PARAMETERS 
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J  -  —5? — F2(z)|!  +  aaeo  (b-a)h1(|,n)  (“ijr)n+1 


12a  L  p  „ 

«  «  -T5?-V*)S.  +  «Sa  h2(l'n>(fr) 


Ac  -  |t-V2(z,|-  +  aeoa  h3(l'n)(5r>n 


«  _  (1-v2 ) PL3  .  PLr  3 (l+\>)  _  3 ( 1-v2 )  _  3v(l-v2) 
An  '  6EI -  +  b~L  — 51 - 132 - 12 


SINGLE  EDGE  0.21(l-v2)P 

CRACKED  PANEL  IN  ~  E 

THREE-POINT  BENDING 

(STB)  A  =  A„  +  A 

n  c 


SUPPORTING  ELASTIC  FUNCTIONS:  F,  Vj^  and  V2  with  z  = 


F(Z)  = 

1.09  -  1.735z  + 

8.2z2  -  14 . 18z 3  + 

14 . 57z " 

Vx(z) 

=  0.76  -  2 . 28z  + 

3 . 87z2  -  2 . 04z  3  + 

0.66/ (1-z ) 2 

N 

CM 

> 

-  (rr^)2  <5.58  - 

19.57z  +  36 . 82z 2 

-  34 . 94z  3  +  12 . 77z " ) 

where 

z  =  a/b 

THEORETICAL  PLASTIC  LIMIT  LOAD:  P^ 

PLANE  STRESS  CONDITION 
P^  =  0.536oq (b-a) 2/L 

PLANE  STRAIN  CONDITION 
p£  =  0.728oq (b-a) 2/L 

SUPPORTING  PLASTIC  FUNCTIONS:  h1#  h2,  h3 

Plane  stress  condition  summarized  in  Table  5a 
Plane  strain  condition  summarized  in  Table  5b 


Figure  8.  Elastic-Plastic  Parameters  for  the  Single-Edge 

Cracked  Panel  in  Three-Point  Bending  (References 
4,  5  and  17) . 
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AND  h  FOR  THE  PLANE  STRESS  SECP  UNDER  THREE-POINT  BENDING 

(REFERENCE  5) . 
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AND  h  FOR  THE  PLANE  STRAIN  SECP  UNDER  THREE-POINT  BENDING 
J  (REFERENCE  5) . 
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GEOMETRY 


ELASTIC-PLASTIC  PARAMETERS 
waeF2 (z)P2 

iTF* 

4a„V0 (z)P 


+  cto  e  a  (ti~— h  (a/b,n)  (L?) 

0  0  D  1  pi 


P_  x  n+1 
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v  0  \  £. )  r  T  n 

-  +  0  e  a  h  (a/b,n)  (P/PMn 


IdE1  ■  — 0“  “3 

4a  V. (z)P 


A  = 


-  +  ae0a  h3(a/b,n)  (P/Pq) n 


A  =  +  /Tae  L  ( -s— - — ) 

n  Eb  o  2bo 


/IP  .  n 


SINGLE  EDGE  CRACKED 

PANEL  IN  TENSION  (SET)  A  =  Ar  +  Ac 


SUPPORTING  ELASTIC  FUNCTIONS:  F,  Vj_  and  V 2  with  z  - 

F  (z)  =  /(2/ttz)  Tan  (ttz/2)  [0.752  +  2.02z  +  0.37(1  -  Sin  (ttz/2  )} 3  ] 
Sec  (irz/2) 

V^z)  =  [  1.46  +  3.42(1  -  Cos  (ttz/2)  }  ]  (Sec  ttz/2)2 

V2(z)  =  [ z/ (1-z ) 2 ] [ 0 . 99  -  z (1-z) (1.3  -  1 • 2z  +  0.7z2)] 

T 

THEORETICAL  PLASTIC  LIMIT  LOAD:  PQ 

PLANE  STRESS  CONDITION 

PT  =  1.072  n  (b-a) 0 
o  0 

PLANE  STRAIN  CONDITION 

PT  =  L.ii  n  (b-a)  a 
°  /3 


with  n  -  2+1  "  (E=a) 

SUPPORTING  PLASTIC  FUNCTIONS:  hj.,  h2 ,  h3 

Plane  stress  condition  summarized  in  Table  6a 
Plane  strain  condition  summarized  in  Table  6b 


Figure  9.  Elastic-Plastic  Parameters  for  Single-Edge 

Cracked  Panel  in  Tension  (References  4,  5  and  17). 
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AND  FOR  THE  PLANE  STRESS  SECP  IN  TENSION 
(REFERENCE  5). 
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AND  h3  FOR  THE  PLANE  STRAIN  SECP  IN  TENSION 
(REFERENCE  5). 
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SECTION  3 

COMPUTER  PROGRAM  DESCRIPTION 


The  purpose  of  this  section  is  to  describe  the  computer 
program  (EST)  listed  in  Appendix  A.  The  EST  computer  program 
was  written  in  FORTRAN  IV  for  interactive  use  on  the  CYBER  175. 
The  first  subsection  describes  its  function  and  the  major  ele¬ 
ments  of  the  program  using  flow  diagrams.  The  second  subsection 
describes  the  program  output  and  various  methods  for  presenting 
the  output  data.  The  third  subsection  presents  the  input  and  how 
one  should  prepare  the  material  properties  data.  The  final  sub¬ 
section  describes  the  format  of  Appendix  A. 

3.1  PROGRAM  FUNCTION 

The  purpose  of  the  EST  computer  program  is  to  calculate  the 
value  of  the  J-Integral  (J) ,  the  crack  mouth  opening  displacement 
(CMOD  =  6),  and  the  load-line  displacement  (A),  as  well  as  their 
elastic-plastic  components.  The  material  is  assumed  to  behave 
according  to  Ramberg-Osgood  type  of  stress-strain  relationship: 


a 

a 


+  a 


(f»n 

o 


(19) 


where  e  and  o  are  reference  strains  and  stresses  related  by 

oo 

eQ  =  oq/e ,  E  is  the  elastic  modulus,  and  a  and  n  are  power 
hardening  constants.  The  program  calculates  the  parameters 
according  to  the  estimating  procedures  defined  within  subsection 
2.2  and  uses  the  equations  and  data  presented  in  Figures  5 
through  9  and  Tables  2  through  6  in  subsection  2.3.  The  values 
of  functions  h^,  h2  and  h^  are  known  for  the  particular  exponent 
n  and  crack  length  a.  In  EST,  the  intermediate  values  of  h^, 
h2  and  h^  are  evaluated  using  a  two-dimensional  linear  interpo¬ 
lation  algorithm  from  the  tabulated  values. 
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One  modification  of  the  theory  was  incorporated  into  the 
EST  computer  program  to  change  the  level  of  the  theoretical 
limit  load  so  that  it  more  closely  approximated  the  observed 
experimental  results.  The  change  involved  utilizing  a  limit 
load  (P  )  calculated  from  the  product  of  a  constant  k  (dependent 
on  material  and  geometry)  and  the  theoretical  limit  load,  i.e. 

P  =  kP  T  (20) 

o  o 

in  all  those  equations  which  utilized  the  limit  load  in  the 
calculation  of  the  elastic-plastic  parameters.  The  value  of  k 
is  an  input  to  the  computer  program  and  would  typically  be  set 
equal  to  1  for  the  first  set  of  runs.  Subsequent  runs  made  with 
this  interactive  program  would  then  be  made  to  bring  the  theore¬ 
tical  load-displacement  behavior  more  in  line  with  the  observed 
behavior. 


3.2  OVERALL  FLOW  DIAGRAM 

Flow  diagrams  of  the  computer  program  are  shown  in 
Figures  10  and  11.  The  call  to  the  subroutine  RDDATA  will 
initiate  the  entering  of  the  test  conditions  interactively. 
Tabulated  values  of  h^,  h2  and  h^  in  Tables  2  through  6  are 
organized  into  separate  subroutines  in  the  program.  Once  the 
geometry  and  stress  condition  are  interactively  chosen,  the  pro¬ 
gram  calls  the  corresponding  subroutine  with  h^,  h 2  and  h^. 

With  the  input  material  constants  and  crack  length 
information,  the  program  makes  a  call  to  subroutine  TBL2 .  Sub¬ 
routine  TBL2  is  a  module  which  returns  the  linearly  interpolated 
values  of  h.^,  h2  and  h^  to  the  program  for  the  given  n  and  a/b 
condition. 


After  the  interpolation  of  the  corresponding  h1# 
h2  and  h3,  the  program  will  make  calls  to  the  calculating  sub¬ 
routines  JKCAL,  LDSCAL  and  CODCAL.  The  subroutine  JKCAL  is 
written  in  such  a  way  that  for  the  selected  geometry,  it  will 

m 

compute  K,  P^,  Yy  and,  elastic  and  plastic  components  of  the 
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/interactively  input 

GEOMETRY  TYPE,  PLANE 
/STRESS  OR  STRAIN, 

/  AND  MATERIAL 

CALL  A  SUBROUTINE  FOR 
TABULATED  VALUES  OF 
hi,  h2,  h3  (TABLES  2-6) 
DEPENDENT  ON  GEOMETRY 
TYPE  AND  PLANE  STRESS 
OR  STRAIN 


CALL  RDDATA  TO  INTER¬ 
ACTIVELY  ENTER  TEST 
CONDITIONS 


PRINT  TEST 
CONDITIONS 


CALL  TBL2  TO  INTERPOLATE 


hj,  h? 


h3  FOR  GIVEN 
MATERIAL  AND  CRACK 
LENGTH 


PRINT  hi,  h*  h 


!•  n2’  n3 


CALL  JKCAL  TO  CALCULATE 
Pj,.  K,  J  AND  ITS 
COMPONENTS 


CALL  LDSCAL  TO 
CALCULATE  A  AND  ITS 
COMPONENTS 


CALL  CODCAL  TO 
CALCULATE  8  AND 
ITS  COMPONENTS 


Figure  10.  Overall  Flow  Diagram  of  the  Computer  Program. 
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Figure  11a.  Detailed  Flow  Diagram  of  the  Computer  Program, 
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Figure  lib.  Detailed  Flow  Diagram  of  the  Computer  Program. 
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Figure  11c.  Detailed  Flow  Diagram  of  the  Computer  Program, 
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J-integral  and  return  these  values  to  the  main  program.  Sub- 

0 

routines  LDSCAL  and  CODCAL  compute  and  return  the  values  of  An, 
Ap,  A®  and,  6e  and  6P  respectively,  to  the  main  program. 

Once  these  calculations  are  completed  for  the 
given  crack  length  and  load,  calculated  values  are  printed.  Then 
the  load  values  and  crack  length  values  are  increased  and  calcu¬ 
lations  are  performed  as  described  in  the  flow  diagrams  in 
Figures  10  and  11.  This  procedure  continues  until  a/b  exceeds 
0.75. 


3.3  PROGRAM  OUTPUT 

The  program  has  been  organized  so  that  the  output  is 
currently  restricted  to  yield  information  on  only  one  geometry 
and  one  material  for  each  run.  The  major  variable  controlling 
output  is  physical  crack  length;  the  other  parameters  are 
described  as  a  function  of  increasing  load  for  each  crack  length. 
Table  7  has  been  prepared  to  provide  a  listing  of  the  output 
parameter  symbols  and  their  description.  Appendix  B  provides 
the  output  for  a  typical  run.  Table  8  describes  some  of  the 
cross-correlations  which  are  possible  based  on  the  EST  program 
output.  Figure  12  was  prepared  to  schematically  describe  how 
several  of  the  possible  cross-correlations  would  appear.  Alter¬ 
nate  presentation  schemes  are  possible  through  a  slight  recon¬ 
figuration  of  the  algorithm  described  in  the  flow  diagrams. 
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CRACK  LENGTH,  a 


STRESS 

INTENSITY 

FACTOR 

K 


Figure  12 


Sample  of  Program  Output  Presented  in  Schematic  Diagrams 


TABLE  7 


PARAMETER  OUTPUT  AND  THEIR  SYMBOLS 


Computer 

Text 

Units 

Symbol 

Symbol 

Description 

(English/Metric) 

A 

a 

Physical  Crack  Length 

in/m 

AE 

ae 

Effective  Crack  Length 

in/m 

ALPHA 

a 

Hardening  Constant 

— 

CMOD 

6 

Crack  Mouth  Opening 
Displacement 

in/m 

E 

E 

Elastic  (or  Young's) 
Modulus 

ksi/Pa 

EC 

Elastic  Crack  (Load- 
line)  Displacement 

in/m 

EN 

c 

Elastic  No-Crack 
(Load-line) Displace¬ 
ment  Component 

in/m 

EPO 

e 

0 

Reference  Strain 

— 

HI 

h1 (a/b,n) 

J-Integral  Plastic 
Function 

— 

H2 

h2 (a/b,n) 

CMOD  Plastic  Function 

— 

H3 

h3 (a/b,n) 

Plastic  Load-line 
Displacement  Function 

— 

J 

J 

J- Integral 

kip-in/N-m 

K 

K  or  SIF 

Stress  Intensity  Factor  ksi/In/Pa/r 

LC 

k 

Limit  Load  Correction 
Factor 

-- 
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TABLE  7  (Continued) 


PARAMETER  OUTPUT  AND  THEIR  SYMBOLS 


Computer 

Symbol 

Text 

Symbol 

Description 

Units 

(English/Metric' 

LDISP 

A 

Load  Point  Displacement 

in/m 

N 

n 

Ramberg-Osgood  Exponent 

— 

PC 

*cP 

Plastic  Crack  (Load  Line) 
Displacement 

in/m 

PN 

Plastic  No-Crack 
(Load  Line) 

in/m 

P/PO 

p/po 

Ratio  of  Load  to  Limit 
Load 

— 

PZC 

— 

Plastic  Zone  Correction 
Factor  (Monotonic=l ; 
Cyclic=2 ) 

SIGO 

ao 

Yield  Strength 

k  s  i  /  P  a 

SPAN 

2L 

Gauge  Length 

in/m 

SQRT(EJ) 

/ej 

Square  Root  of  (E.J) 

ksi/Tn/Pa/in 

SQRT ( J/ JE) 

/j/Jti 

Square  Root  of  (J/Je) 

— 

TH 

Thickness  of  the 

Specimen 

in/m 

# 
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TABLE  3 

POTENTIAL  CROSS-CORRELATIONS  POSSIBLE  WITH  PROGRAM  OUTPUT 


Primary 

Independent 

Variable 

(Abcissa) 


Load-line 

Displacement 

•  Total 

•  Elastic 

•  Plastic 

Crack  Length 


Crack  length 


Load-Line 

Displacement 


Dependent 

Variable 

(Ordinate) 


Secondary 

Independent 

Variable 

(Multiple  Curves) 


BASIC  EXPERIMENTAL  PARAMETERS 


Load 


Crack  Length 


Component 

Component 


Elastic  Compliance  - 

Maximum  Load-line 

Displacement  Load 

Load  to  Limit  Load 

Ratio  Load 

Maximum  Crack  Mouth 

Opening  Displacement  Load 


FATIGUE  AND  FRACTURE  RELATED  PARAMETERS 


Plastic  Zone  Size 

Effective  Crack  Length 

Stress  Intensity  Factor 
J-Integral  (HRR  Field) 

•  /jTe 

•  smz 

Crack  Tip  Opening 
Displacement 


Load,  Field 
Parameters 

Load,  Field 
Parameters 

Load 

Load,  Displacement, 

Load  to  Limit  Load 
Ratio 


Load 


J-Integral  (HRR  Field) 

•  /J.E 

•  73^ 


Crack  Length 
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3.4  PROGRAM  INPUT  REQUIREMENTS 


Necessary  inputs  to  the  EST  computer  program  are  material, 
structural  geometry,  and  loading  properties.  A  summary  listing 
of  the  input  parameters  and  their  symbols  are  presented  in 
Table  9  and  a  sample  input  listing  is  provided  in  Appendix  C. 

As  can  be  noted  from  this  listing,  the  material  hardening  para¬ 
meter  a  does  not  appear.  The  reason  for  this  is  that  EST 
calculates  a  (ALPHA)  internally  assuming  (1)  that  the  constitutive 
equation  is  the  Ramberg-Osgood  stress-strain  relationship: 

=  -2-  +  a  (— ) n  (21) 

£  0  0 
0  0  0 

(2)  that  the  reference  stresses  and  strains  are  related  by 


=  Ee 


(3)  that 


aeo  =  «(-#>  = 


0.002 


(22) 


and  (4)  that  the  hardening  exponent  n  is  obtained  by  least  squares 
fitting  procedures.  Equation  22  was  derived  with  the  use  of 
Figure  13  and  an  analysis  of  Equation  21  for  the  condition 

o=o  which  shows 

o 


ey=e|  =£0(1+a) 
J  0=0' 


(23) 


When  the  EST  program  is  compiled  and  run,  the  program  prompts 
the  user  for  necessary  information.  Input  of  the  geometry,  as 
well  as  selection  of  plane  stress  or  plane  strain  condition  allow 
the  program  to  select  the  appropriate  table  of  values  for  h^,  hj 


and  h^  functions  presented  in  subsection  2.3.  The  EST  program 
calculates  the  EPFM  parameters  for  different  loads,  starting  from 
the  user  selected  minimum  load  and  ending  at  the  user  selected 


maximum  load  with  load  increments  selected  by  the  user  for  a 
given  initial  a/b  ratio.  By  a  suitable  selection  of  crack  length 
increment,  it  is  possible  to  calculate  the  EPFM  parameters  of 
different  crack  lengths  starting  from  the  initial  (a/b)  for 


different  loads. 
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TABLE  9 

EST  COMPUTER  PROGRAM  INPUTS 


EST 

SYMBOL 

TEXT 

SYMBOL 

PROPERTY 

DEFAULT 

VALUE 

ENGLISH/ 

METRIC 

UNITS 

SIGO 

0 

0 

Material  0.2  Percent  Yield 
Strength 

ksi/Pa 

N 

n 

Ramberg-Osgood  Exponent 

— 

E 

E 

Elastic  Modulus 

ksi/Pa 

LC 

k 

Correction  for  Limit  Load 

1.0 

— 

PZC 

Plastic  Zone  Correction 

Factor 

1.0 

—  —  — 

— 

Plane  Stress/Plane  Strain 
Condition 

— 

— 

Specimen  Geometry  (CCP,  CT, 
DCP ,  STB,  SET) 

— 

w 

2b 

Structural  Width 

in/m 

— 

L 

Half  Span 

in/m 

TH 

Thickness 

in/m 

A/B 

a/b 

Initial  Crack  Length/ 

Width  Ratio 

0.25 

— 

Crack  Length  Increment 

in/m 

Load/Displacement 

PMAX 

Maximum  Load 

kips/N 

PMIN 

Minimum  Load 

0.0 

kips/N 

Load  Increment 

kips/N 

DISMAX 

Maximum  Load-Line 

Displacement 

in/m 
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Figure  13.  Material  Stress-Strain  Curve  Describing  the 
Values  of  a  and  eo* 
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In  the  input  prompts,  the  plastic  zone  correction  factor 
(PC)  is  taken  as  either  1  or  2.  The  value  PC  is  used  in  Equation 
16  as  a  multiplication  factor  on  the  yield  strength  oq;  when 
PC  =  1  the  loading  is  monotonically  icreasing,  when  PC  =  2  the 
loading  is  cyclic. 

There  are  several  restrictions  or  bounds  on  the  calculations 
that  result  from  input  or  default  conditions.  The  most  obvious 
bounds  are  the  limits  on  the  load  or  load-line  displacement  values 
that  are  created  by  input.  However,  because  the  tabular  values 
of  the  hlf  h2,  and  h3  functions  are  limited  to  the  range  of  a/b 
between  0.25  and  0.75,  the  EST  program  frustrates  any  attempt 
to  compute  the  value  of  these  functions  outside  the  allowable 
a/b  range  by  resetting  upper  and  lower  limits.  Another  default 
condition  occurs  whenever  the  handling  exponent  (n)  exceeds  20 
in  order  to  ensure  that  the  h^,  h2,  h^  functions  are  not  extra¬ 
polated  beyond  the  tabulated  values  given  in  subsection  2.3. 
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SECTION  4 


SUMMARY 

A  review  of  the  theoretical  basis  for  the  estimation  scheme 

12  3  6 

proposed  by  Hutchinson  and  co-workers  '  and  Shih  and  co-workers  ' 
was  presented  in  previous  sections.  A  computer  program  was  written 
to  implement  the  above  estimation  scheme.  Certain  modifications 
were  incorporated  into  the  program  such  that  the  limit  load 
behavior  computed  by  the  program  is  compatible  with  the  experi¬ 
mentally  observed  limit  load  behavior.  The  computer  program  can 
be  used  to  calculate  the  elastic-plastic  behavior  of  five  different 
specimen  geometries.  Possible  useful  outputs  from  the  program 
are  given  in  Table  8. 


* 
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APPENDIX  A 


A  LISTING  OF  THE  FORTRAN  IV 
PROGRAM  EST 

DESIGNED  TO  RUN  INTERACTIVELY 
ON  THE  CDC  CYBER  175 
COMPUTER  SYSTEM 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

11 


2100 


100 

2200 


2500 

2600 

200 


PROGRAM  EST( INPUT, OUTPUT, TAPE6 ,TAPE1=INPUT) 

COMMON/ A/E , EA , PLCOR , SIG0 , XLL , B , FACTOR , XN , ALPHA , EP0 , TH , PP 1 , I PL , 

$  A,H1,H2,H3, PP , ZZ , ZZE , P0 , AE ,XKSQ , XJ , BETA , PI 
COMMON/ B/DELTE , DELTP , DELTA 

COMMON/C/DELCE , DELCP , DELNE , DELNP , DELC , DELN , DELT , DELE , DELP 
COMMON/D/X ( 10 )  , Y ( 10 )  ,Z1(100)  ,Z2(100)  ,Z3(100)  , I FL AG , METENG 
COMMON/E/  PMIN , PMAX , PIN , ZZ0 ,AIN,DISMAX 
COMMON/F/NUMAB ,NUMN , I TYPE 
COMMON/X/PPPO ,XJJE 

PROGRAM  EST  CALCULATES  J,  DISPLACEMENT,  AND  DELTA 
FOR  THE  FOLLOWING  SPECIMENS: 

1)  CENTER  CRACKED  PANEL  (CCP)  -  PLANE  STRAIN  OR  PLANE  STRESS 

2)  COMPACT  TENSION  SPECIMEN  (CT)  PLANE  STRAIN 

3)  DOUBLE- EDGE  CRACK  PANEL  (DCP)  PLANE  STRAIN 

4)  SINGLE-EDGE  CRACK  PANEL  IN  3 -POINT  BENDING  -  PLANE  STRAIN 

5)  SINGLE-EDGE  CRACK  PANEL  IN  TENSION  (SET)  -  PLANE  STRAIN 

THE  VALUES  OF  DISPLACEMENT  VS.  J,  LOAD,  AND  DELTA  ARE 
WRITEN  TO  A  FILE  USED  TO  GENERATE  PLOTS  (UNIT  7) . 


REWIND  6 
XNU=0 .3 
IFLAG=0 
PRINT2100 

FORMAT (*0ENTER  SPECIMEN  GEOMETRY*,/, 

$  *  (CCP ,CT , DCP , STB , SET) == |  *) 

READ (1, 100 ) JTYPE 
FORMAT (A3) 

PRINT2200 

FORMAT (*  PLANE  STRESS  (0)  OR  PLANE  STRAIN 
READ* , IPL 
PRINT2500 
FORMAT (*  MATERIAL 
READ (1,2600) MAT 
FORMAT (A10) 

WRITE (6, 200) MAT 
FORMAT (*  MATERIAL 
IF ( JTYPE. EQ.3HCCP. AND. IPL. EQ.0) 

$  CALL  CCPSTS 

IF ( JTYPE. EQ.3HCCP. AND. IPL. EQ.l) 

$  CALL  CCPSTN 
IF ( JTYPE. EQ.3HCT 
$  CALL  CTSTS 
IF ( JTYPE. EQ.3HCT 
$  CALL  CTSTN 
IF (JTYPE .EQ . 3HDCP 
$  CALL  DCPSTS 
IF ( JTYPE. EQ.3HDCP 
$  CALL  DCPSTN 
IF ( JTYPE. EQ.3HSTB  .AND. 

$  CALL  STBSTS 
IF ( JTYPE. EQ.3HSTB 
$  CALL  STBSTN 


(1)  ?  *) 


==l  *) 


* , A10 ) 


.AND. 


.AND. 


IPL.EQ.0) 
IPL. EQ.l) 


.AND.  IPL.EQ.0) 


.AND. 


IPL. EQ.l) 
IPL.EQ.0) 


.AND.  IPL. EQ.l) 
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IF ( JTYPE . EQ . 3HSET  .AND.  IPL.EQ.O) 

$  CALL  SETSTS 

IF ( JTYPE. EQ.3HSET  .AND.  IPL.EQ.l) 

$  CALL  SETSTN 
IF ( IFLAG .EQ . 0 ) PRINT2300 
2300  FORMAT (*0 INVALID  GEOMETRY*) 

IF ( IFLAG. EQ.O)  GO  TO  11 
C 

C  READ  VARIABLES  INTERACTIVELY 

C 

CALL  RDDATA 

A0=ZZ0*B 

A=A0 

BETA=2 . 0 

IF ( IPL  .NE.  0)  BETA=6 . 0 
WRITE (6 ,300) 

300  FORMAT (6X, *TH* , 9X, *SPAN* , 11X, *N* , 13X , *SIG0* , 7X , *PZC* , 14X, 

$  *E* , 8X , * WIDTH* ,10X,*LC*,9X, * ALPHA* , 10X , *EP0* ,/) 

EA=E 

IF ( IPL . NE . 0 ) EA=E/ ( 1 . -XNU**2 ) 

EP0=SIG0/E 
ALPHA=0 . 002/EP0 
SPAN=2*XLL 
B2=B 

I F ( I TYPE . EQ . 1 . OR . I TYPE . EQ . 3 )  B2=B*2 

WRI TE ( 6 , 4  0  0 ) TH , SPAN , XN , SIGO , PLCOR , E , B2 , FACTOR , ALPHA , EPO 
400  FORMAT (10G13. 4) 

12  ZZ=A/B 

IF (XN .GT .20.0)  XN=20 . 0 
WRITE (6,1300) XN , ZZ 

1300  FORMAT ( *0  X  =  * ,G14 . 4 , 5X , *A/B  =  *,G14.4) 

C 

C  INTERPOLATE  FOR  VALUES  OF  H1,H2,H3 

C 

CALL  TBL2 (X , Y , Z1 ,NUMN , NUMAB , NUMN , XN , ZZ ,H1 , IER) 

CALL  TBL2 (X ,Y , Z2 , NUMN , NUMAB , NUMN ,XN , ZZ , H2 , IER) 

CALL  TBL2 ( X , Y , Z3 , NUMN , NUMAB , NUMN , XN , ZZ , H3 , IER) 

WRITE (6 ,1400) H1,H2,H3 

1400  FORMAT ( *  HI  =  * ,F12 . 4 , 3X , *H2  =  * , F12 . 4 , 3X , *H3  =  *,F12.4) 

IF (METENG  .EQ.  1HM)  WRITE (6 , 1500) 

IF (METENG  .EQ.  1HE)  WRITE (6 , 1600 ) 

1500  FORMAT (*0ALL  RESULTS  IN  METRIC  UNITS*,/) 

1600  FORMAT (*0ALL  RESULTS  IN  ENGLISH  UNITS*,/) 

IF (ITYPE.EQ. 1) WRITE (6,500) 

500  FORMAT (T4,*2A*,T8,*2A/W*,T15,*2AE*,T21,*P/P0*,T31, 

$  *LOAD* ,T42,*K*,T53,*J*,T59, *SQRT (EJ) * ,T68 , *SQRT ( J/JE) * , 

$  T8 1 , *CMOD* ,T9 1 , *LDISP* ,T10 2 , *EN* ,T111 , *PN* ,T119 , *EC* , T128 , *PC* ) 
IF(ITYPE.NE.l)  WRITE (6,600) 

600  FORMAT (T5,*A* , T9 , *A/W* ,T15 , *AE* ,T21 , *P/P0  * ,T31 , 

$  *LOAD* ,T42,*K*,T53,*J*,T59, *SQRT (EJ) * ,T68 , *SQRT ( J/JE) * , 

$  T81 , *CMOD* ,T91 , *LDISP* ,T102 , *EN* ,T111 , *PN* ,T119 , *EC* ,T128 , *PC*) 
PP1=PMIN 

13  CONTINUE 
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CALCULATE  J  AND  K 
CALL  JKCAL 

CALCULATE  LOAD-POINT  DISPLACEMENT 
CALL  LDSCAL 

CALCULATE  CRACK-OPENING  DISPLACEMENT 
CALL  CODCAL 
C 

XK=SQRT ( XKSQ ) 

XXJ=SQRT(XJ*EA) 

XJJE=SQRT (XJJE) 

AE2=AE 

A2=A 

IF ( ITYPE  .EQ.  1)  A2=A*2 
IF ( ITYPE  .EQ.  1)  AE2=AE*2 

WRITE ( 6 , 7  00 ) A2 , Z  Z , AE2 , PPPO , PP1 , XK , XJ , XXJ , XJJE , DELTA , DELT , DELNE , 
$DELNP , DELCE , DELCP 

700  FORMAT ( 1H  , 2F5 . 3 , F6 . 3 , E10 . 3 , F9 . 3 , 3E10 . 3 , F10 . 7 , E10 . 3 / E9 . 3 , E10 . 3 , 
$  3E9.3) 

IF (DELT  .GE.  DISMAX)  GO  TO  15 
PP1=PP1+PIN 

IF  (PP1  .LT.  PMAX)  GO  TO  13 
15  CONTINUE 
A=A+AIN 

IF  (A/B  .LE.  3  ,/4 . )  GO  TO  12 
STOP 
END 
C 
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SUBROUTINE  RDDATA 

**  INITIALIZES  VARIABLES  AND  PROMPTS  FOR  USER  INPUT  ** 

COMMON/A/E , EA , PLCOR , SIGO , XLL , B , FACTOR , XN , ALPHA , EPO ,TH,PP1,IPL, 
$A ,  HI , H2 ,H3,PP,ZZ,ZZE,P0 ,  AE , XKSQ ,XJ , BETA , PI 
COMMON/D/X ( 10 ) ,Y(10) ,Z1(100) ,Z2(100) , Z3 ( 100 ) , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 
COMMON/E/  PMIN , PMAX , PIN , ZZO , AIN , DISMAX 
DATA  PMIN/0 . 0/ , PIN/ . 25/ , PMAX/5 . 3/ , ZZO/ . 25/ ,XN/36 ./ , 

$  SIGO/33E3/,PLCOR/l./,E/17E6/,XLL/2./,B/1.995/, 

$  TH/0 .498//AIN/0 . 1/ , FACTOR/ 1./ ,DISMAX/0 .01/ 

PRINT1600 

READ (1,2000) METENG 

PRINT100 ,  PMIN 

READ (1 ,2000) KANS 

IF(KANS.NE.IHY) READ*, PMIN 

PRINT200 , PIN 

READ (1,2000) KANS 

IF(KANS.NE.IHY) READ*, PIN 

PRINT300 ,PMAX 

READ (1,2000) KANS 

I F ( KANS . NE . 1HY ) READ* , PMAX 

PRINT500 , ZZO 

READ (1,2000) KANS 

IF(KANS.NE.IHY) READ*, ZZO 

PRINT600 ,XN 

READ (1,2000) KANS 

I F ( K ANS . NE . 1HY ) READ  * , XN 

PRINT700 ,SIG0 

READ (1,2000) KANS 

IF ( KANS . NE . 1HY) READ* , SIGO 

PRINT800 , PLCOR 

READ (1,2000) KANS 

I F ( KANS . NE . 1HY ) READ* , PLCOR 

PRINT900 ,E 

READ (1,2000) KANS 

IF(KANS.NE.IHY) READ* ,E 

PRINT1000 , XLL 

READ (1,2000) KANS 

I F ( KAN  S . NE . 1HY ) READ  * , XLL 

PRINT1100 , B 

READ (1,2000) KANS 

IF (KANS .NE . 1HY) READ* ,B 

IF ( I TYPE  .EQ.  1  .OR.  I TYPE  .EQ.  3)B=B/2.0 

PRINT1200 ,TH 

READ (1,2000) KANS 

IF(KANS.NE.IHY) READ* ,TH 

PRINT1300 , AIN 

READ (1,2000) KANS 

IF (KANS.NE. 1HY) READ*, AIN 

PRINT1400 , FACTOR 

READ (1,2000) KANS 

IF ( KANS . NE . 1HY) READ* , FACTOR 

PRINT1500 , DISMAX 

READ (1,2000) KANS 

IF (KANS.NE. 1HY) READ*, DISMAX 
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100  FORMAT ( *0IF  THE  FOLLOWING  VALUES  ARE  CORRECT,*,/, 

$*  RESPOND  WITH  A  "Y" .  IF  THEY  ARE  INCORRECT,*,/, 

$*  ANSWER  "N",  THEN  ENTER  THE  CORRECT  VALUE*,/, 

$*0MINIMUM  LOAD  (KIPS/N)  =  *,G10.5,*  (Y/N) ?  *) 

200  FORMAT (*  LOAD  INCREMENT  (KIPS/N)  =  *,G11.5,*  (Y/N)?  *) 

300  FORMAT ( *  MAXIMUM  LOAD(KIPS/N)  =  *,G7.2,*  (Y/N)?  *) 

500  FORMAT ( *  INITIAL  A/B  =  *,G11.5,*  (Y/N)?  *) 

600  FORMAT (*  RAMBERG -OSGOOD  EXPONENT  -  N  =  *,G11.5,*  (Y/N)?  *) 

700  FORMAT (*  YIELD  STRENGTH  -  SIG0  (KSI/PA)  =  *,G14.5,*  (Y/N)?  *) 

800  FORMAT (*  PLASTIC  ZONE  CORRECTION  FACTOR  -  PZC=  *,G14.5,*  (Y/N)?  *) 
900  FORMAT ( *  ELASTIC  MODULUS  -  E  (KSI/PA)  =  *,G14.5,*  (Y/N)?  *) 

1000  FORMAT ( *  HALF  SPAN  (IN/M)  =  *,G11.5,*  (Y/N)?  *) 

1100  FORMAT ( *  WIDTH  -  W  (IN/M)  =  *,G11.5,*  (Y/N)?  *) 

1200  FORMAT (*  THICKNESS  -  TH  (IN/M)  =  *,G12.5,*  (Y/N)?  *) 

1300  FORMAT (*  CRACK  LENGTH  INCREMENT  (IN/M)  =  *,G10.5,*  (Y/N)?  *) 

1400  FORMAT(*  CORR.  FACTOR  FOR  LIMIT  LOAD  -  LC  =  *,G10.5,*  (Y/N)?  *) 
1500  FORMAT(*  MAXIMUM  DISPLACEMENT  (IN/M)  =  *,G10.5,*  (Y/N)?  *) 

1600  FORMAT (*  E(NGLISH)  OR  M(ETRIC)  ?  *) 

2000  FORMAT (Al) 

RETURN 

END 


54 


oooonooooooooooooooonooooooooooooooonooooooonoooo 


SUBROUTINE  TBL2  (X,  Y,  Z,  NX,  NY,  NDX ,  XO,  YO,  ZO,  IER) 

TBL2  FROM  THE  COMPUTER  CENTER  LIBRARY  OF  6600  ROUTINES 
SUBROUTINE  TBL2 

PERFORMS  DOUBLE  LINEAR  INTERPOLATION  FROM  A  TABLE  OF  VALUES 

OF  Z  VERSES  X  AND  Y. 

USAGE 

CALL  TBL2  (X , Y , Z ,NX ,NDX , XO , YO , ZO , IER) 

DESCRIPTION  OF  PARAMETERS 

X  VALUES  OF  THE  FIRST  INDEPENDENT  VARIABLE  IN  INCREASING 
ORDER 

Y  VALUES  OF  THE  SECOND  INDEPENDENT  VARIABLE  IN  INCREASING 
ORDER 

Z  VALUES  OF  THE  DEPENDENT  VARIABLE.  Z(I,J)  IS  THE  VALUE 
OF  Z  CORRESPONDING  TO  X(I),  Y(J)  AND  IS  REFERENCED  AS 
Z( (J-l) *NDX+I) 

NX  NUMBER  OF  ENTRIES  IN  THE  X  ARRAY 

NY  NUMBER  OF  ENTRIES  IN  THE  Y  ARRAY 

NDX  DIMENSION  IN  THE  X  DIRECTION  OF  THE  Z  ARRAY  (NDX .GE .NX) 

XO  X  COORDINATE  OF  POINT  AT  WHICH  INTERPOLATION  IS  DESIRED 

YO  Y  COORDINATE  OF  POINT  AT  WHICH  INTERPOLATION  IS  DESIRED 

ZO  THE  COMPUTED  INTERPOLATED  VALUE 

IER  *  0  INTERPOLATION  SUCESSFULLY  PERFORMED 

=  1  EXTRAPOLATION  SUCESSFULLY  PERFORMED 

=  2  ERROR  CONDITION.  EITHER  TWO  ADJACENT  INDEPENDENT 
VARIABLES  ARE  NOT  IN  INCREASING  ORDER  OR  NX  IS  GREATER 
THAN  NDX 


METHOD 

XO  AND  YO  ARE  LOCATED  RELATIVE  TO  THE  X  AND  Y  ARRAYS.  IF 
EXTRAPOLATION  IS  NOT  NECESSARY,  REPEATED  BISECTION  OF  THE 
INDICES  OF  THE  LISTS  OF  X  AND  Y  VALUES  IS  CARRIED  OUT  UNTIL 
BOTH  ARGUMENTS  ARE  ISOLATED  BETWEEN  PAIRS  OF  CONSECUTIVE 
VALUES.  AFTER  LOCATION  OF  THE  NEAREST  PAIR  OF  ARGUMENTS  FOR 
EACH  OF  THE  TWO  INDEPENDENT  VARIABLES,  DOUBLE  LINEAR  INTER¬ 
POLATION  (EXTRAPOLATION)  IS  PERFORMED. 


EXAMPLE  -  GIVEN  THE  FUNCTION  Z 
AS  A  ONE  DIMENSIONAL  OR  TWO 


FOLLOWS  Z(l) 

=  Z  ( 1 , 1) 

2 

1 

Z(2) 

»  Z  ( 2 , 1) 

s 

2 

WHERE 

THE  DIMENSION  OF 

z 

IS 

OR 

Z  (1) 

=  Z(l,l) 

= 

1 

z  (2) 

=  Z  ( 2 , 1) 

= 

2 

WHERE 

THE  DIMENSION  OF 

z 

IS 

=  X+2Y ,  X=1 , 2  AND  Y=0 , I  STORAGE 
DIMENSIONAL  ARRAY  COULD  BE  AS 
Z  (3)  =  Z  (1,2)  =  3 
Z  ( 4)  =  Z  (  2 , 2)  =  4 
Z ( 4 )  OR  Z (2 ,2)  AND  NDX  =  2 
Z (5)  =  Z(l,2)  =  3 
Z  (6 )  =  Z  ( 2 , 2 )  *  4 
Z ( 16 )  OR  Z  ( 4 , 4 )  AND  NDX  =  4. 


C 


DIMENSION  X(l),  Y  ( 1)  ,  Z(l) 
IF  (NX.GT.NDX)  GO  TO  12 
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C  DETERMINE  WHETHER  EXTRAPOLATION  IS  NECESSARY  IN  THE  X  DIRECTION 

11  =  1 

IF  (XO.LT.X(l) )  GO  TO  4 
IF  (XO.GT.X(NX) )  GO  TO  3 

LOCATE  XO  WITHIN  X ( I )  ARRAY  FOR  INTERPOLATION 
IER  =  0 

12  =  NX 

1  IF  ( (12-11) .LT.2)  GO  TO  5 

I  =  ( I 1+12) /2 

IF  (XO.LT.X (I) )  GO  TO  2 

II  =  I 
GO  TO  1 

2  12  =  I 
GO  TO  1 

SET  INDICES  AND  FLAG  FOR  EXTRAPOLATION 

3  II  =  NX-1 

4  IER  =  1 

DETERMINE  WHETHER  EXTRAPOLATION  IS  NECESSARY  IN  THE  Y  DIRECTION 

5  J1  =  1 

IF  (YO.LT.Y(l) )  GO  TO  9 
IF  (YO.GT.Y(NY) )  GO  TO  8 

LOCATE  YO  WITHIN  Y(I)  ARRAY  FOR  INTERPOLATION 
J2  =  NY 

6  IF  ( (J2-J1) .LT.2)  GO  TO  10 
J  =  ( J1+J2) /2 
IF  (YO.LT.Y(J) )  GO  TO  7 
J1  =  J 
GO  TO  6 

7  J2  =  J 
GO  TO  6 

SET  INDICES  AND  FLAG  FOR  EXTRAPOLATION 

8  J1  =  NY-1 

9  IER  =  1 

COMPUTE  ZO  USING  BIVARIATE  INTERPOLATION 

10  DIV  =  (X(I1+1)  -  X (II) )  *  ( Y ( Jl+1)  -Y(J1)) 

IF  (DIV)  12,12, 11 

11  Ill  =  ( Jl-1) *NDX  +  II 
112  =  Ill  +  NDX 
X20  =  X(I1+1)  -  XO 
X01  =  XO  -  X(I1) 

ZO  =  (Y2-YO) ( (X2-XO) Zll  +  (XO-X1) Z21)  + 

(YO-Y1) ( (X2-XO) Z12  +  (XO-X1) Z22)  /  (X2-X1) (Y2-Y1) 

ZO  =  ( (Y (Jl+1) -YO) *(X20*Z(I11)  +  X01*Z(I11+1) )  + 

+  (YO-Y(Jl)  ) *(X20*Z(I12)  +  X01*Z(I12+1) ) )  /DIV 

RETURN 

12  IER  =  2 
RETURN 
END 
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non 


SUBROUTINE  JKCAL 

**  CALCULATES  ELASTIC  AND  PLASTIC  COMPONENTS  OF  J  ** 

COMMON/ A/E , EA , PLCOR , SIGO , XLL , B , FACTOR , XN , ALPHA , EPO , TH , PP1 , I PL , 

$  A , HI , H2 , H3 , PP , ZZ , Z ZE , PO , AE , XKSQ , XJ , BETA , PI 
COMMON/X/PPPO ,  XJJE 
COMMON/F/NUMAB , NUMN , I TYPE 
FSQCCP(Z)  =  l./(COS( (3.14159/2.) *Z) ) 

FSQCT(Z)  =*( ( (2+Z) *(0.886+4. 64*Z-13 . 32*Z**2+14 . 72*Z**3 
$  -5.6*Z**4) ) **2) / ( Z* ( 1-Z) **3) 

FSQSTB (Z)=(1.09-1.735*Z+8.2*Z**2-14.18*Z**3+14.57*Z**4) **2 
FSQDCP (Z)=( 1.122-0. 561*Z-0. 205*Z**2+0. 471*Z**3-0. 19*Z**4) **2/ (1-Z) 
FSQST(Z) =(2/(PI*Z) ) *TAN(PI*Z/2.) * ( (0 . 752+2 . 02*Z+0 . 37 
$  * (1-SIN (PI *Z/2. ) ) **3) /COS (PI*Z/2 . ) )**2 
PI=3. 14159 
PP=PP1/TH 
ZZ=A/B 

GO  TO  (1,2,3,4,5)ITYPE 
C  CCP  SPECIMEN 

1  P0=2.*SIG0*(B-A) ‘FACTOR 

IF ( IPL  .NE.  0)  P0=2.31*SIG0*(B-A) *FACTOR 
XKSQ= (PP**2*PI*A*FSQCCP ( ZZ) ) / ( 4 . *B**2) 

PHI=1 ./ ( 1. + ( PP/PO ) **2) 

RY= (XN-1 . ) *XKSQ/ ( BETA*PI * ( XN+1 . ) * ( SIGO*PLCOR) **2) 

AE=A+PHI*RY 

ZZE-AE/B 

XJE= (PI*AE*FSQCCP ( ZZE) * (PP**2)  ) 

$  / (EA*4 . *B**2) 

XJP=ALPHA*SIGO*EPO*A* ( (B-A) /B) *H1* (PP/PO ) ** (XN+1 . ) 

GO  TO  6 

C  CT  SPECIMEN 

2  XNETA=SQRT ( ( ( 2 . * A) / ( B-A) ) **2+ ( ( 4 . *A) / (B-A)  ) 

$  +2  . )  -  (  ( 2  .  *A)  /  (B-A)  +1. ) 

P0=1.071*XNETA* (B-A) *SIGO*FACTOR 

IF ( IPL  .NE.  0)  P0=1.455*XNETA* (B-A) *SIG0*FACTOR 

XKSQ=(PP**2*A*FSQCT(ZZ) )/(B**2) 

PHI-1./ (1.+ (PP/PO ) **2) 

RY= (XN-1. ) *XKSQ/(BETA*PI* (XN+1.) * ( SIG0*PLCOR) **2) 

AE=A+PHI*RY 

ZZE=AE/B 

XJE= (AE*FSQCT ( ZZE) * (PP**2) ) / (EA*B**2) 

XJP=ALPHA*SIG0*EP0* (B-A) *H1* (PP/PO) ** (XN+1) 

GO  TO  6 

C  DCP  SPECIMEN 

3  P0=2 . 31*SIG0* (B-A) *FACTOR 
IF(IPL.NE.O)  P0=5.94*SIG0* (B-A) * FACTOR 
XKSQ= (PP**2*PI*A*FSQDCP ( ZZ) )/(2.*B) **2 
PHI=1./ (1.+ (PP/PO ) **2) 

RY= (XN-1. ) *XKSQ/ ( BETA*PI * (XN+1. ) * ( SIGO*PLCOR) **2) 

AE=A+PHI*RY 

ZZE=AE/B 

XJE= (PI *AE*FSQDCP ( ZZE) * (PP**2) )/(EA*B**2) 

XJP=ALPHA*SIGO*EPO* (B-A) *H1* (PP/PO )** (XN+1 . ) 

GO  TO  6 
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c  STB  SPECIMEN 

4  P0= (0 . 536*SIG0* (B-A) **2/XLL) * FACTOR 
IF(IPL.NE.O)  P0=(0.728*SIG0* (B-A) **2/XLL) * FACTOR 
XKSQ=9.*PP**2*XLL**2*PI*A*FSQSTB(ZZ)/( (B) **4) 

PHI-1. /(l.+(PP/P0) **2 ) 

RY= (XN-1. ) *XKSQ/ (BETA*PI* (XN+1. ) * (SIGO*PLCOR) **2) 

AE=A+PHI*RY 

ZZE=AE/B 

XJE= (9 . *PI *AE*XLL**2*FSQSTB ( ZZE) *PP**2)/(EA* (B) **4) 
XJP=ALPHA*SIGO*EPO* (B-A) *H1* (PP/PO) ** (XN+1. ) 

GO  TO  6 

C  SET  SPECIMEN 

5  XNETA=SQRT ( 1+ ( A/ (B-A) ) **2) - (A/ (B-A) ) 

P0=1.455*XNETA* (B-A) *SIGO*FACTOR 

IF ( IPL . EQ . 0 ) P0  =  1 . 07  2*XNETA* ( B-A) *SIG0  *FACTOR 
XKSQ=(PP**2*PI*A*FSQST(ZZ) )/(B**2) 

PHI-1.0/ (1.0+ (PP/PO) **2) 

RY= (XN-1 . 0 ) *XKSQ/(BETA*PI* (XN+1.0) * (SIGO*PLCOR) **2) 

AE=A+PHI *RY 

ZZE=AE/B 

XJE= (PI*AE*FSQST (ZZE) * (PP**2) )/(EA*B**2) 
XJP=ALPHA*SIGO*EPO*A* ( (B-A)/B) *H1* (PP/PO) ** (XN+1. ) 

6  XJ=XJE+XJP 
XJEE= (XKSQ/EA) 

IF(XJEE.EQ.O.O)  XJJE-0.0 

IF ( XJEE . NE .0.0)  XJJE= (XJ/XJEE) 

PPPO-PP/PO 

RETURN 

END 
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OO  cn  o  **  O  LJ  n  too  M  o  ooooooo 


SUBROUTINE  LDSCAL 

COMMON/ A/E , EA , PLCOR ,  SIGO , XLL , B , FACTOR , XN , ALPHA , EPO , TH , PP1 , IPL , 
$  A , HI , H2 , H3 ,PP,ZZ,ZZE,PO , AE , XKSQ , XJ , BETA , PI 
COMMON/C/DELCE , DELCP , DELNE , DELNP , DELC , DELN , DELT , DELE , DELP 
COMMON/F/NUMAB , NUMN , I TYPE 

**  LOAD  POINT  DISP  CALC.  FOR  ALL  SPECIMENS  ** 

DELCE  =  ELASTIC  CRACKED  COMPONENT 
DELCP  =  PLASTIC  CRACKED  COMPONENT 
DELNE  =  ELASTIC  NO  CRACKED  COMPONENT 
DELNP  =  PLASTIC  NO  CRACKED  COMPONENT 

V2CCP (2) =-1.071+0 .25*Z-0.357*Z**2+0.121*Z**3 
$  -0.047*Z**4+0.008*Z**5-1.071*(1./Z)  MALOG(l.-Z)  ) 

V2CT( Z)®(11.062*Z**3-27.209*Z**2+27 . 977*Z+0 . 995) / 

$  ((1-Z)**2) 

V2DCP (Z)=(2/(PI*Z) ) *(0.0629-0. 061*COS(PI*Z/2.) **4 
$  -0 . 00 19*COS (PI*Z/2 . ) **8+ALOG ( 1/ (COS ( PI *Z/2 . ) ) ) ) 

V2STB (Z)=(Z/(1-Z) ) **2*(5. 58-19. 57*Z+36.82*Z**2 
$  -34.94*Z**3+12.77*Z**4) 

V2ST ( Z)  =  ( Z/ (1-Z) **2) *(0.99-Z*(l-Z) *( 1.3-1. 2*Z+0.7*Z**2)  ) 

XNU=0 . 3 

DELCP=ALPHA*EP0*A*H3* (PP/P0) **XN 
GO  TO  (1,2, 3, 4, 5) ITYPE 
CCP  SPECIMEN 

DELCE = ( 2 . *AE*V2CCP ( ZZE) *PP) / (B*EA) 

DELNE= ( PP*XLL) / ( B*E) 

DELNP=1 . 7  321*ALPHA*EP0*XLL* ( ( 0 . 4 33*PP) / (B*SIG0 ) ) 

$  **XN 

GO  TO  6 

CT  SPECIMEN 
DELCE=V2CT (ZZE) *PP/EA 
DELNE-0 .0 
DELNP=0 . 0 
GO  TO  6 

DCP  SPECIMEN 

DELCE=2*AE*V2DCP ( ZZE) *PP/(B*EA) 

DELCP=ALPHA*EP0* (B-A) *H3* (PP/P0) **XN 
DELNE® ( PP*XLL) / (B*E) 

DELNP® 1 . 7321*ALPHA*EP0*XLL* ( (0.433*PP)/ (B*SIG0 ) ) 

GO  TO  6 

STB  SPECIMEN 

DELCE® (6*XLL**2*V2STB( ZZE) *PP) / (EA*B**2) 

DELCP® ALPHA*EP0 *A*H3* (PP/P0 ) **XN 

DELNE® ( (l.-XNU**2)*PP*XLL**3)/(E*0.5*B**3)+( (PP*XLL) / (B*E) ) 

#  * (1.5* (l.+XNU) -0 . 3* ( 1 . -XNU**2 ) -0 . 75*XNU* ( 1 . -XNU**2) ) 

*  -(0.21*(1.-XNU**2) *PP)/E 
DELNP®0 . 

GO  TO  6 

SET  SPECIMEN 

DELCE® ( 4 . *AE*V2ST ( ZZE) *PP) / (EA*B) 

DELCP=ALPHA*EP0*A*H3* (PP/P0) **XN 
DELNE* ( 2 . *PP*XLL) / (B*E) 

DELNP*1 . 7321*ALPHA*EP0*XLL* ( (0 . 866*PP) / ( 2 . *B*SIG0 ) ) **XN 
DELC  *  TOTAL  DISPLACEMENT  DUE  TO  CRACK 
DELN  =  TOTAL  DISPLACEMENT  DUE  TO  NO  CRACK 
6  DELC*DELCE+DELCP 
DELN=DELNE+DELNP 
DELT*DELC+DELN 
DELE*DELNE+DELCE 
DELP*DELNP+DELCP 
RETURN  59 

END 


UOUUU  CJ  CJ  (N  CJ  cj  ** 


SUBROUTINE  CODCAL 


**  CALCULATES  CRACK-OPENING  DISPLACEMENT  FOR  ALL  SPECIMENS  ** 
DELTE  =  ELASTIC  COMPONENT 
DELTP  =  PLASTIC  COMPONENT 

COMMON/A/E , EA , PLCOR , SIGO , XLL , B , FACTOR , XN , ALPHA , EPO ,TH,PP1,IPL, 

$  A , HI , H2 , H3 , PP , ZZ , ZZE , PO , AE , XKSQ , XJ , BETA , PI 
COMMON/B/DELTE , DELTP , DELTA 

COMMON/C/DELCE , DELCP , DELNE  f  DELNP , DELC , DELN , DELT , DELE , DELP 
COMMON/F/NUMAB , NUMN , I TYPE 

VICCP(Z) =-0. 071-0. 535*Z+0.169*Z**2+0.02*Z**3 
$  -1.071* (l./Z) * (ALOG ( 1 . -Z) ) 

V1CT(Z) =(57.694*Z**3-83.166*Z**2+43.315*Z+5.435)/ 

$( (1-Z) **2) 

VIDCP(Z) =(2/(PI*Z) ) * ( (0 .459*SIN(PI*Z/2.) ) - ( 0 . 06 5* ( SIN (PI *Z/2 . ) ) 
$  **3) -(0.007* (SIN (PI*Z/2. ) ) **5 ) + ( ALOG ( ( 1/ (COS ( PI *Z/2 . ) ) ) 

$  +SQRT ( ( 1/ (COS ( PI *Z/2 . ) ) ) **2-1.)  )  )  ) 

VISTB(Z) =0 .76-2.28*Z+3 .87*Z**2-2.04*Z**3+0 .66/ (1-Z) **2 
VIST ( Z) = (1. 46+3. 42* (1-COS (PI*Z/2.) ) ) /COS (PI*Z/2 . ) **2 
GO  TO  (1,2,3,4,5)ITYPE 
CCP  SPECIMEN 

1  DELTE= (2. *AE*V1CCP(ZZE) *PP)/(B*EA) 

DELTP=ALPHA*EP0*A*H2* (PP/P0) **XN 
GO  TO  6 

CT  SPECIMEN 
DELTE=V1CT ( ZZE) *PP/EA 
DELTP=ALPHA*EP0*A*H2* (PP/P0) **XN 
GO  TO  6 

DCP  SPECIMEN 

DELTE=2 . *AE*V1DCP (ZZE) *PP/(B*EA) 

DELTP=ALPHA*EP0  * (B-A) *H2* (PP/P0) **XN 
GO  TO  6 

STB  SPECIMEN 

DELTE=12 . * AE*XLL*V1STB ( ZZE) *PP/ ( B**2*EA) 

DELTP= ALPHA* EPO *A*H2* (PP/P0) **XN 
GO  TO  6 
C  SET  SPECIMEN 

5  DELTE=4 . *AE*V1ST (ZZE) *PP/(B*EA) 

DELTP=ALPHA*EP0  * A*H2* (PP/P0) **XN 

6  DELTA=DELTE+DELTP 
RETURN 

END 
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SUBROUTINE  CCPSTS 

COMMON/D/X ( 10 )  f Y ( 10)  ,Z1(100)  ,Z2(100)  ,Z3(100)  , I FLAG , METENG 
COMMON/F/NUMAB , NUMN  , I TYPE 

DATA  X/l. , 1.5,2. ,3. ,5. ,7. ,10. ,13. ,16. ,20./ 


DATA  Y/0.125 

,0.2,0 

.25,0. 

375,0. 

5,0.625 

i,0.75 

,3*0.0/ 

DATA  Zl/ 

$  2.80,  3.23, 

3.54, 

4.00, 

4.52, 

4.76, 

4.86, 

3.25, 

3.95, 

3.65 

$  2.63,  3.00, 

3.15, 

3.42, 

3.62, 

3.65, 

3.43, 

2.95, 

2.89, 

2.55 

$  2.54,  2.82, 

2.97, 

3.14, 

3.19, 

3.11, 

2.90, 

2.65, 

2.47, 

2.20 

$  2.34,  2.40, 

2.53, 

2.52, 

2.35, 

2.17, 

1.95, 

1.77, 

1.61, 

1.43 

$  2.21,  2.24, 

2.20, 

2.06, 

1.81, 

1.63, 

1.43, 

1.30, 

1.17, 

1.00 

$  2.11,  2.07, 

1.91, 

1.69, 

1.41, 

1.22, 

1.01, 

0.85, 

0.71, 

0.57 

$  2.07,  1.89, 

1.71, 

1.46, 

1.21, 

1.08, 

0.96, 

0.75, 

0.65, 

0.53 

$  30*0.0/ 

DATA  Z2/ 

$  3.53,  3.86, 

4.11, 

4.64, 

4.83, 

4.94, 

4.89, 

4.40, 

3.37, 

3.11 

$  3.28,  3.45, 

3.60, 

3.73, 

3.61, 

3.55, 

3.25, 

2.95, 

2.57, 

2.48 

$  3.12,  3.24, 

3.29, 

3.30, 

3.15, 

2.93, 

2.59, 

2.29, 

2.08, 

1.81 

$  2.71,  2.65, 

2.62, 

2.41, 

2.03, 

1.75, 

1.47, 

1.28, 

1.13, 

0.99 

$  2.34,  2.18, 

2.01, 

1.70, 

1.30, 

1.07, 

1.87, 

0.76, 

0.67, 

0.56 

$  1.97,  1.72, 

1.46, 

1.13, 

0.79, 

0.62, 

0.47, 

0.38, 

0.31, 

0.26 

$  1.61,  1.22, 

0.97, 

0.69, 

0.45, 

0.36, 

0.29, 

0.22, 

0.18, 

0.15 

$  30*0.0/ 

DATA  Z3/ 

$  0.35,  0.49, 

0.64, 

0.95, 

1.54, 

2.05, 

2.63, 

3.35, 

3.38, 

3.11 

$  0.52,  0.69, 

0.90, 

1.24, 

1.77, 

2.09, 

2.37, 

2.55, 

2.57, 

2.48 

$  0.61,  0.82, 

1.01, 

1.35, 

1.83, 

2.08, 

2.19, 

2.12, 

2.01, 

1.79 

$  0.81,  1.00, 

1.19, 

1.43, 

1.59, 

1.57, 

1.43, 

1.27, 

1.13, 

0.99 

$  0.93,  1.09, 

1.19, 

1.26, 

1.18, 

1.04, 

0.87, 

0.76, 

0.67, 

0.56 

$  0.98,  1.07, 

1.05, 

0.97, 

0.76, 

0.62, 

0.48, 

0.39, 

0.32, 

0.27 

$  0.93,  0.89, 
$  30*0.0/ 

0.80, 

0.64, 

0.45, 

0.36, 

0.29, 

0.22, 

0.18, 

i 

0.15 

WRITE(6 ,100) 

100  FORMAT (*  CENTER  CRACK  PANEL  -  PLANE  STRESS*,//) 

ITYPE=1 

NUMAB=7 

NUMN=10 

I FLAG =1 

RETURN 

END 
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SUBROUTINE  CCPSTN 

COMMON/D/X ( 10 ) , Y ( 10) ,Z1(100) ,Z2(100) ,Z3(100) , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX ( 7 ) ,XY(7) ,XZ1(49) ,XZ2(49) ,XZ3(49) 

DATA  XX/  1., 1*5, 2. ,3. ,5. ,7. ,10./ 


DATA  XY/  0.25,0 
DATA  XZ1/ 

.3,0.4, 

0.5, 0.6 

o 

• 

** 

o 

• 

7  5/ 

$ 

2.570, 

2.900, 

3.110, 

3.350, 

3.490, 

3.430, 

3.230, 

$ 

2.450, 

2.720, 

2.890, 

3.100, 

3.110, 

2.950, 

2.720, 

$ 

2.300, 

2.450, 

2.550, 

2.600, 

2.450, 

2.230, 

1.950, 

$ 

2.190, 

2.250, 

2.270, 

2.180, 

1.930, 

1.710, 

1.440, 

$ 

2.120, 

2.060, 

2.030, 

1.830, 

1.530, 

1.350, 

1.060, 

$ 

2.100, 

1.900, 

1.870, 

1.630, 

1.230, 

1.110, 

0.840, 

$  2.100,  1.850, 
DATA  XZ2/ 

1.800, 

1.570, 

1.240, 

1.040  , 

0.799/ 

$ 

2.790, 

2.990, 

3.090, 

3.140, 

3.000, 

2.790, 

2.470, 

$ 

2.630, 

2.800, 

2.820, 

2.790, 

2.560, 

2.280, 

1.970, 

$ 

2.350, 

2.380, 

2.320, 

2.160, 

1.800, 

1.520, 

1.270, 

$ 

2.090, 

1.990, 

1.870, 

1.610, 

1.230, 

0.996, 

0.776, 

$ 

1.810, 

1.580  , 

1.420, 

1.140, 

1.780, 

0.610, 

0.450, 

$ 

1.540, 

1.230, 

1.050, 

0.630, 

0.490, 

0.360, 

0.250, 

$  1.400,  1.080, 
DATA  XZ3/ 

0.899, 

0.637, 

0.401, 

0.295, 

0.204/ 

$ 

0.548, 

0.752, 

0.942, 

1.270, 

1.730, 

1.970, 

2.070, 

$ 

0.600, 

0.830, 

1.010, 

1.290, 

1.620, 

1.770, 

1.750, 

$ 

0.710, 

0.920, 

1.070, 

1.260, 

1.380, 

1.350, 

1.190, 

$ 

0.798, 

0.949, 

1.060, 

1.150, 

1.100, 

0.959, 

0.771, 

$ 

0.810, 

0.900, 

1.000, 

0.970, 

0.810, 

0.600, 

0.600, 

$ 

0.820, 

0.820, 

0.780, 

0.720, 

0.530, 

0.360, 

0.360, 

$ 

0.814, 

0.763, 

0.733, 

0.593, 

0.399, 

0,294, 

0.204/ 

WRITE (6,100) 

100  FORMAT (*  CENTER  CRACK  PANEL  -  PLANE  STRAIN*,//) 

DO  10  1=1,7 
X(I) =XX(I) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,49 
Z1(I) =XZ1(I) 

Z2 (I) =XZ2 (I) 

Z3 ( I ) =XZ3 ( I ) 

20  CONTINUE 

ITYPE=1 
IFLAG=1 
NUMAB=7 
NUMN=7 
RETURN 
END 
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SUBROUTINE  CTSTS 

COMMON/D/X (10) , Y ( 10 ) ,Z1(100) ,Z2(100) ,Z3(100) , I FLAG , METENG 


COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX(10) ,XY(10) ,XZ1(72) ,XZ2(72) ,XZ3(72) 
DATA  XX/1. ,2. ,3. ,5. ,7. ,10. ,13. ,16  .  ,20.  ,0.0/ 


DATA  XY/ 0.25 

,0.30,0 

.375,0 

• 

o 

in 

rf 

• 

• 

o 

o 

in 

525,0.' 

75,1.00 

o 

• 

o 

to 

DATA  XZ1/ 

$  1.61,  1.46, 

1.28, 

1.06, 

0.90, 

0.73, 

0.60, 

0.51, 

0.39, 

$  1.60,  1.35, 

1.24, 

0.94, 

0.79, 

0.59, 

0.49, 

0.39, 

0.28, 

$  1.55,  1.25, 

1.05, 

0.80, 

0.65, 

0.48, 

0.38, 

0.28, 

0.22, 

$  1.47,  1.11, 

0.93, 

0.78, 

0.58, 

0.46, 

0.35, 

0.24, 

0.20, 

$  1.40,  1.08, 

0.90, 

0.69, 

0.56, 

0.44, 

0.36, 

0.30, 

0.24, 

$  1.27,  1.03, 

0.88, 

0.69, 

0.59, 

0.49, 

0.42, 

0.37, 

0.31, 

$  1.23,  0.98, 

0.83, 

0.68, 

0.60, 

0.51, 

0.43, 

0.37, 

0.31, 

$  1.13,  1.01, 

0.78, 

0.68, 

0.65, 

0.62, 

0.49, 

0.47, 

0.42/ 

DATA  XZ2/ 

$17.55,12.04, 

10.71, 

8.74, 

7.32, 

5.74, 

4.63, 

3.75, 

2.92, 

$15.10,10.27, 

8.70, 

6.80, 

5.10, 

3.90, 

3.31, 

2.20, 

1.71, 

$12.41,  8.20, 

6.54, 

4.56, 

3.45, 

2.44, 

1.83, 

1.36, 

1.02, 

$10.31,  6.50, 

5.02, 

3.32, 

2.51, 

1.73, 

1.32, 

1.09, 

0.82, 

$  9.16,  5.67, 

4.21, 

2.80, 

2.12, 

1.57, 

1.25, 

1.03, 

0.81, 

$  7.47,  4.48, 

3.35, 

2.37, 

1.92, 

1.54, 

1.29, 

1.12, 

0.93, 

$  6.25,  3.78, 

2.89, 

2.14, 

1.78, 

1.44, 

1.20, 

1.03, 

0.86  , 

$  5.29,  3.54, 

2.41, 

1.91, 

1.73, 

1.59, 

1.23, 

1.17, 

1.03/ 

DATA  XZ3/ 

$  9.67,  8.00, 

7.21, 

5.94, 

5.00, 

3.95, 

3.19, 

2.59, 

2.02, 

$  8.80,  7.10, 

6.10, 

6.88, 

3.60, 

2.82, 

2.22, 

1.69, 

1.30, 

$  7.80,  5.73, 

4.62, 

3.25, 

2.48, 

1.77, 

1.33, 

0.99, 

0.75, 

$  6.68,  4.73, 

3.68, 

2.42, 

1.81, 

1.27, 

1.00, 

1.80, 

0.60, 

$  6.29,  4.15, 

3.11, 

2.09, 

1.59, 

1.18, 

0.94, 

0.77, 

0.61, 

$  5.42,  3.38, 

2.54, 

1.80, 

1.47, 

1.18, 

0.99, 

0.85, 

0.71, 

$  4.77,  2.92, 

2.24, 

1.66, 

1.38, 

1.12, 

0.94, 

0.80, 

0.67, 

$  4.23,  2.83, 

1.93, 

1.52, 

1.39, 

1.27, 

0.99, 

0.93, 

0.82/ 

WRITE(6 ,100) 

100  FORMAT (*  COMPACT  TENSION  SPECIMEN - PLANE  STRESS*,//) 

ITYPE=2 
NUMAB=8 
NUMN=9 
DO  1  1-1,10 
X(I)-XX(I) 

1  Y ( I )  =XY ( I ) 

DO  2  1=1,72 
Z1(I)=XZ1(I) 

Z2 ( I ) =XZ2 ( I ) 

2  Z3 ( I ) =XZ3 ( I ) 

I FLAG =1 
RETURN 

END 
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10 


20 


100 


SUBROUTINE  CTSTN 

COMMON/D/X (10)  , Y (10) ,Z1(100) ,Z2(100) ,Z3(100) , I FLAG , METENG 
COMMON/ F/NUMAB , NUMN , I TYPE 

DIMENSION  XX(10) ,XY(10) ,XZ1(63) ,XZ2(63) ,XZ3(63) 


DATA  XX/1. ,2. ,3 

•  ,  5  •  ,  7  . 

,10. ,13 

. ,16. , 

20.  ,0.0/ 

3*0.0/ 

DATA  XY/0.25,0. 
DATA  XZ1/ 

3,0.375 

,0.45,0 

.5,0.625,0.75, 

1.325 

$  2.227, 

2.048, 

1.783, 

1.475, 

1.334, 

1.248, 

1.253, 

$  2.200, 

1.870, 

1.540, 

1.190, 

0.960, 

0.780, 

0.700, 

0.500 

$  2.148, 

1.716 , 

1.392, 

0.970, 

0.693, 

0.443, 

0.276, 

0.176 

$  2.040, 

1.580, 

1.280, 

0.900, 

0.540, 

0.390, 

0.280, 

0.200 

$  1.935, 

1.509, 

1.242, 

0.919, 

0.685, 

0.461, 

0.314, 

0.216 

$  1.763, 

1.449, 

1.237, 

0.974, 

0.752, 

0.602, 

0.459, 

0.347 

$  1.709, 

1.424, 

1.263, 

1.033, 

0.864, 

0.717, 

0.575, 

0.448 

1.566, 

0.520, 

0.098, 

0.060, 

0.132, 

0.248, 

0.345/ 


DATA  XZ2/ 

$17.883,12.481,11.675,10.788,10.538,10.745,11.460,12.570,14.563 


$15,850,10.430  , 

8.900 

$12,644,  8.176, 

6.521 

$10,550,  6.650, 

5.030 

$  9.327,  5.846, 

4.304 

$  7.612,  4.572, 

3.423 

$  6.370,  3.948, 

3.179 

DATA  XZ3/ 

$  9.852,  8.506, 

8.170 

$  9.000,  7.100, 

4.170 

$  7.944,  5.760, 

4.643 

$  7.020,  4.730, 

3.560 

$  6.406,  4.268, 

3.157 

$  5.521,  3.431, 

3.179 

$  4.857,  3.048, 

2.456 

I FLAG =1 

ITYPE=2 

DO  10  1=1,10 

X(I)=XX(I) 

Y ( I ) =XY ( I ) 

DO  20  1=1,63 

Z1(I)=XZ1(I) 

Z2 ( I ) =XZ2 ( I ) 

Z  3 ( I ) =XZ3 ( I ) 

NUMAB=7 

6.020, 

5.250, 

3.840 

4.319, 

2.970, 

1.744 

3.180, 

2.200, 

0.820 

2.747, 

1.912, 

1.199 

2.359, 

1.810, 

1.319 

2.337, 

1.876, 

1.441 

7.774, 

7.706, 

7.942 

4.660, 

3.540, 

3.050 

3.103, 

2.139, 

1.292 

2.360, 

1.630, 

0.910 

2.024, 

1.413, 

0.888 

1.787, 

1.373, 

1.000 

1.807, 

1.450, 

1.114 

3.700, 

7.750, 

9.200, 

1.102, 

0.686, 

0.370, 

0.830, 

0.520, 

0.310, 

0.788, 

0.530  , 

0.317, 

0.983, 

0.749, 

0.485, 

1.124, 

0.887, 

0.685/ 

8.517, 

9.371, 

10.887, 

2.600, 

5.700, 

6.700, 

0.793, 

0.494, 

0.266, 

0.580, 

0.360, 

0.200, 

0.585, 

0.393, 

0.236, 

0.746, 

0.568, 

0.368, 

0.869, 

0.686, 

0.514/ 

//) 


NUMN=9 
WRITE (6 ,100) 

FORMAT (*  COMPACT  TENSION  -  PLANE  STRAIN* 


RETURN 

END 
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SUBROUTINE  DCPSTN 

COMMON/ D/X ( 10 ) ,  Y  ( 10) ,Z1(100) ,Z2(100) ,23(100) , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX (9) ,XY(9) ,XZ1(45) ,XZ2(45) ,XZ3(45) 

DATA  XX/1.,  2.,  3.,  5.,  7.,  10.,  13.,  16.,  20./ 

DATA  XY/0.25,  0.375,  0.5,  0.625,  0.75,  4*0.0/ 

DATA  XZ1/ 

$  5.01,128.4,28.67,155.43,749.07,7530.8,75618. ,741530. ,156. E5, 

$  6.41,14.56 ,30.42,120.99,461.24,3354.2,20  374. ,168810. ,2278. E3, 
$  7.31,14.59,27.07,87.53,275.57,1505.9,8109.4,43102. ,343780. , 

$  7.87,13.40,21.62,53.51,128.26,455.26,1558.4,5376.8,28037. , 

$  8.19,12.4,17.72,33.06,55.12,105.32,221.79,518. ,894.59/ 

DATA  XZ2/ 

$  3.33,7.41,16.55,78.46,362.04,3507.8,34044. ,326360. ,6719300. , 

$  5.00,9.61,18.64,67.83,244.2,1676.3,11491. ,80300. ,1062000. , 

$  6.78,11.22,18.94,53.76,155.16,784.89,4008.9,20719. ,173640. , 

$  8.79,11.91,16.73,34.22,73.4,238.5,783.6,2623.4,13141. , 

$  11.37,12.21,14.07,20.72,31.92,58.76,117.99,262.33,473.82/ 
DATA  XZ3/ 

$  0.57,2.2,6.82,49.33,297.75,3720.9,43022. ,459780 ., 10421 .E3 , 

$  1.27,4.22,11.17,59.32,266.63,2211.5,16698. ,123400. ,1693100. , 

$  2.28,6.39,14.29,55.49,187.17,1051.1,5610.1,29420. ,244640. , 

$  3.73,8.28,14.83,38.63,91.15,310.61,1030.7,3458.7,17316. , 

$  5.91,9.89,13.89,23.82,38.01,70.98,142.31,315.28,575.54/ 

DO  10  1=1,9 
X ( I ) =XX ( I ) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,45 
Z1(I)=XZ1(I) 

Z2 ( I ) =XZ2 ( I ) 

Z3 ( I ) =XZ3 ( I ) 

20  CONTINUE 

WRITE(6 ,100) 

100  FORMAT (*  DOUBLE-EDGE  CRACK  PANEL  ---  PLANE  STRAIN*,//) 

I TYPE =3 

IFLAG=1 

NUMAB=5 

NUMN=9 

RETURN 

END 
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SUBROUTINE  STBSTN 

COMMON/D/X ( 10 )  ,  Y  ( 10)  ,Z1(100)  ,22(100)  ,23  (100)  , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX (10) ,XY(10) ,X21(63) ,XZ2(63) ,X23(63) 

DATA  XX/1. ,2. ,3. ,5. ,7. ,10. ,13. ,16. ,20. ,0.0/ 

DATA  XY/0. 2 5, 0.3, 0.35, 0.4, 0.5, 0.6 ,0.75,3*0.0/ 


DATA  X21/ 

$  1.195,  1.034, 

0.930, 

0.765, 

0.633, 

0.523, 

0.396, 

0.303,  0.215, 

$ 

1.270,  1.040, 

0.920, 

0.740, 

0.570, 

0.450, 

0.320, 

0.240,  0.150, 

$ 

1.340,  1.050, 

0.910, 

0.720, 

0.540, 

0.410, 

0.280, 

0.200,  0.120, 

$ 

1.370,  1.070, 

0.910, 

0.700, 

0.510, 

0.370, 

0.240, 

0.160,  0.090, 

$ 

1.398,  1.094, 

0.922, 

0.675, 

0.495, 

0.331, 

0.211, 

0.135,  0.074, 

$ 

1.360,  1.100, 

0.920, 

0.670, 

0.480, 

0.320, 

0.250, 

0.110,  0.150, 

$ 

1.208,  1.145, 

0.974, 

0.693, 

0.500  , 

0.348, 

0.223, 

0.140,  0.075/ 

DATA  X22/ 

$  5.799,  4.665, 

4.006, 

3.080, 

2.454, 

1.934, 

1.446, 

1.088,  0.758, 

$ 

5.520,  4.180, 

3.450, 

2.670, 

2.060, 

1.580, 

1.100, 

0.830,  0.550, 

$ 

5.300,  4.620, 

3.120, 

2.330, 

1.760, 

1.300, 

0.870, 

0.650,  0.410, 

$ 

5.140,  3.630, 

2.860, 

2.050, 

1.520, 

1.080, 

0.700, 

0.520,  0.300, 

$ 

4.869,  3.283, 

2.527, 

1.686, 

1.192, 

0.773, 

0.480, 

0.304,  0.165, 

$ 

4.650,  3.030, 

2.300, 

1.440, 

0.980, 

0.630, 

0.380, 

0.250,  0.110, 

$ 

4.474,  2.754, 

2.096, 

1.361, 

0.936, 

0.618, 

0.388, 

0.239,  0.128/ 

DATA  X23/ 

$  4.083,10.099, 

8.413, 

5.864, 

4.466, 

3.421, 

2.542, 

1.901,  1.318, 

$ 

4.280,  7.600, 

6.350, 

4.290, 

3.340, 

2.470, 

1.590, 

1.250,  0.900, 

$ 

4.440,  6.330, 

5.140, 

3.540, 

2.740, 

1.920, 

1.230, 

0.930,  0.650, 

$ 

4.560,  5.500, 

4.420, 

3.050, 

2.290, 

1.550, 

0.980, 

0.700,  0.470, 

$ 

4.687,  4.442, 

3.509, 

2.353, 

1.663, 

1.079, 

0.669, 

0.424,  0.230, 

$ 

4.670,  3.800, 

2.300, 

1.900, 

1.320, 

1.830, 

0.510, 

0.280,  0.140, 

$ 

1.491,  3.159, 

2.407, 

1.557, 

1.068, 

0.704, 

0.441, 

0.272,  0.144/ 

DO  10  1=1,10 
X(I)=XX(I) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,63 
21 (I) =X21 (I) 

22(1) =X22 ( I ) 

23(1) =X23 ( I ) 

20  CONTINUE 

WRITE (6 ,100) 

100  FORMAT (*  SINGLE  EDGE  CRACK  PANEL  IN  3-POINT  BENDING  - * 

$  *  PLANE  STRAIN*,//) 

ITYPE=4 

IFLAG=1 

NUMAB=7 

NUMN =9 

RETURN 

END 
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SUBROUTINE  SETSTN 

COMMON/ A/E , EA , PLCOR , SIGO , XLL , B , FACTOR , XN , ALPHA , EPO ,TH,PP1,IPL, 
$A,H1,H2,H3,PP,ZZ,ZZE,P0,AE,XKSQ,XJ,BETA,PI 
COMMON/D/X (10) ,Y(10) ,Z1(100) ,Z2(100) ,Z3(100) , I FLAG , METENG 
COMMON/E/ PMI N  , PMAX , PIN , ZZO , AIN , DISMAX 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX ( 7 ) ,XY(7) ,XZ1(42) ,XZ2(42) ,XZ3(42) 


DATA  XX/1. ,2. , 

3. ,5. ,7 

• 

O 

• 

o 

.0/ 

DATA  XY/0.25,0 

.3,0.35 

,0.4,0. 

5, 0.6,0 

.75/ 

DATA  XZ1/ 

$  4.338, 

4.768, 

4.639, 

3.815, 

3.056, 

2.170, 

$  4.200, 

3.880, 

3.250, 

2.520, 

2.033, 

1.490, 

$  4.020, 

3.340, 

2.650, 

1.900, 

1.450, 

1.010, 

$  3.820, 

3.030, 

2.230, 

1.460, 

1.040, 

0.660, 

$  3.398, 

2.302, 

1.694, 

0.928, 

0.514, 

0.213, 

$  2.960, 

1.910, 

1.400, 

1.690, 

0.330, 

0.110, 

$  2.342, 

1.607, 

1.245, 

0.769, 

0.477, 

0.233/ 

DATA  XZ2/ 

$  4.756, 

4.559, 

4.281, 

3.391, 

2.639, 

1.808, 

$  4.660, 

3.670, 

3.350, 

2.310, 

1.790, 

1.260, 

$  4.590, 

3.280, 

2.790, 

1.750, 

1.270, 

0.870, 

$  4.530, 

3.030, 

2.230, 

1.870, 

0.910, 

0.580, 

$  4.447, 

2.765, 

1.888, 

0.954, 

0.507, 

0.204, 

$  4.880, 

2.520, 

1.400, 

0.810, 

0.380, 

0.150, 

$  4.316, 

2.515, 

1.789, 

1.027, 

0.619, 

0.296/ 

DATA  XZ3/ 

$10,270, 

7.635, 

5.874, 

3.695, 

2.483, 

1.496, 

$  8.150, 

4.120, 

3.000, 

1.930, 

1.560, 

0.900, 

$  6.050, 

3.000, 

1.880, 

1.290, 

1.030, 

0.580, 

$  4.600, 

2.370, 

1.350, 

0.880, 

0.590, 

0.350, 

$  3.151, 

1.537, 

.9125, 

.4172, 

.2151, 

.0851, 

$  2.520, 

1.150, 

0.750, 

0.500, 

0.060, 

0.040, 

$  2.018, 

1.105, 

.7655, 

.4349, 

.2617, 

.1254/ 

WRITE (6 ,100) 

100  FORMAT (*  SINGLE  EDGE  CRACK  PANEL  IN  TENSION  -  * 

$  * PLANE  STRAIN*,//) 

ITYPE=5 
I FLAG =1 
NUMAB=7 
NUMN=6 
DO  10  1*1,7 
X  (I ) *XX ( I ) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1*1,42 
Z1(I)=XZ1(I) 

Z2 ( I ) =XZ2 ( I ) 

Z3 ( I ) =XZ3 ( I ) 

20  CONTINUE 

RETURN 
END 
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SUBROUTINE  DCPSTS 

COMMON/D/X ( 10 ) ,Y(10) ,Z1(100) ,Z2(100) ,Z3(100) , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX ( 9 ) ,XY(9) ,XZ1(45) ,XZ2(45) ,XZ3(45) 

DATA  XX/I. ,2. ,3. ,5. ,7. ,10. ,13. ,16. » 20./ 


DATA  XY/ 0.25,0. 

37,0.5 

,0.62,0. 

75,0.0 

o 

• 

o 

** 

o 

• 

o 

,0.0/ 

DATA  XZ1/ 

$  1.011,  1.226, 

1.356, 

1.483, 

1.543, 

1.578, 

1.594, 

1.591, 

1.588, 

$  1.293,  1.418, 

1.427, 

1.341, 

1.237, 

1.094, 

0.970, 

0.873, 

0.674, 

$  1.475,  1.466, 

1.378, 

1.168, 

1.010, 

0.845, 

0.732, 

0.625, 

0.208, 

$  1.586,  1.454, 

1.284, 

1.038, 

0.882, 

0.737, 

0.649, 

0.466, 

0.020, 

$  1.652,  1.425, 

1.118, 

0.979, 

0.833, 

0.701, 

0.630, 

0.297, 

0.000/ 

DATA  XZ2/ 

$  1.726,  1.819, 

1.886, 

1.917, 

1.905, 

1.853, 

1.802, 

1.746, 

1.700, 

$  2.594,  2.393, 

2.221, 

1.864, 

1.588, 

1.283, 

1.068, 

0.922, 

0.709, 

$  3.514,  2.821, 

2.337, 

1.670, 

1.277, 

0.944, 

0.762, 

0.630, 

0.232, 

$  4.559,  3.145, 

2.318, 

1.449, 

1.061, 

0.790, 

0.657, 

0.473, 

0.028, 

$  5.896,  3.371, 

2.214, 

1.297, 

0.966, 

0.741, 

0.634, 

0.312, 

0.000/ 

DATA  XZ3/ 

$  0.295,  0.537, 

0.770, 

1.169, 

1.490, 

1.815, 

2.022, 

2.124, 

2.198, 

$  0.658,  1.037, 

1.295, 

1.520, 

1.547, 

1.412, 

1.227, 

1.068, 

0.829, 

$  1.184,  1.581, 

1.691, 

1.563, 

1.320, 

0.080, 

0.808, 

0.662, 

0.265, 

$  1.932,  2.138, 

1.950, 

1.444, 

1.094, 

0.809, 

0.664, 

0.487, 

0.032, 

$  3.063,  2.670, 

2.061, 

1.314, 

0.978, 

0.747, 

0.638, 

0.318, 

0.000/ 

WRITE (6 ,100) 

FORMAT ( *  DOUBLE  EDGE 

CRACK  PANEL  -- 

-  PLAIN 

STRESS 

*,//) 

ITYPE=3 

IFLAG=1 
NUMAB=5 
NUMN=9 
DO  10  1=1,9 
X(I)=XX(I) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,45 
Z1(I) =XZ1(I) 
Z2 (I) =XZ2 (I) 
Z3 (I)=XZ3 (I) 
20  CONTINUE 

RETURN 
END 
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SUBROUTINE  STBSTS 

COMMON/D/X ( 10 )  ,Y(10)  ,Z1(100)  ,Z2(100)  ,Z3(100)  , I FLAG , METENG 
COMMON/F/NUMAB , NUMN  , I TYPE 

DIMENSION  XX(9) ,  XY (9) ,XZ1 (45) ,XZ2 (45) ,XZ3 (45) 

DATA  XX/1. ,2. ,3. ,5. ,7. ,10. ,13. ,16. ,20./ 


DATA  XY/0.25,0 

.37,0.5 

t  0  •  6  2  1 0 

.75,0.0 

o 

• 

o 

o 

• 

o 

,0.0/ 

DATA  XZ1/ 

$  0.869,  0.731, 

0.629, 

0.479, 

0.370, 

0.246, 

0.174, 

0.117, 

0.059, 

$  0.963,  0.797, 

0.680, 

0.527, 

0.418, 

0.307, 

0.232, 

0.174, 

0.105, 

$  1.019,  0.767, 

0.621, 

0.453, 

0.324, 

0.202, 

0.128, 

0.081, 

0.030, 

$  1.051,  0.786, 

0.649, 

0.494, 

0.357, 

0.235, 

0.173, 

0.105, 

0.047, 

$  1.067,  0.786, 

0.643, 

0.474, 

0.343, 

0.230, 

0.167, 

0.110, 

0.044/ 

DATA  XZ2/ 

$  5.690,  4.503, 

3.680, 

2.614, 

1.947, 

1.290, 

0.897, 

0.603, 

0.307, 

$  5.085,  3.732, 

2.929, 

2.071, 

1.580, 

1.134, 

0.841, 

0.626, 

0.381, 

$  4.768,  3.120, 

2.320, 

1.547, 

1.077, 

0.655, 

0.410, 

0.259, 

0.097, 

$  4.551,  2.830, 

2.118, 

1.455, 

1.023, 

0.656, 

0.472, 

0.286, 

0.130, 

$  4.385,  2.656, 

1.967, 

1.329, 

0.928, 

0.601, 

0.427, 

0.280, 

0.114/ 

DATA  XZ3/ 

$  4.007,  8.812, 

7.189, 

4.731, 

3.388, 

2.204, 

1.517, 

1.012, 

0.508, 

$  4.420,  5.533, 

4.482, 

3.172, 

2.409, 

1.726, 

1.277, 

0.948, 

0.575, 

$  4.604,  4.085, 

3.092, 

2.081, 

1.442, 

0.874, 

0.545, 

0.344, 

0.129, 

$  4.617,  3.434, 

2.599, 

1.794, 

1.258, 

0.803, 

0.577, 

0.349, 

0.158, 

$  4.394,  3.012, 

2.235, 

1.510, 

1.052, 

0.680, 

0.483, 

0.316, 

0.129/ 

DO  10  1=1,9 
X(I)=XX(I) 

Y ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,45 
Z1(I) =XZ1(I) 

Z2 (I) =XZ2 (I) 

Z3 (I) =XZ3 (I) 

20  CONTINUE 

WRITE (6 ,100) 

100  FORMAT (*  SINGE  EDGE  CRACK  PANEL  IN  3-POINT  BENDING  ---* 
$  * PLANE  STRESS*,//) 

ITYPE=4 

I FLAG =1 

NUMAB=5 

NUMN=9 

RETURN 

END 
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SUBROUTINE  SETSTS 

COMMON/D/X ( 10 )  ,  Y (10)  ,Z1(100)  ,Z2(100)  ,Z3(100)  , I FLAG , METENG 
COMMON/F/NUMAB , NUMN , I TYPE 

DIMENSION  XX ( 9 ) ,XY(9) ,XZ1(45) ,XZ2(45) ,XZ3(45) 

DATA  XX/I. ,2. ,3. ,5. ,7. ,10. ,13. ,16 . ,20./ 


DATA  XY/0.25,0 

.37,0.5 

,0.62,0 

.75,0.0 

• 

o 

w 

o 

• 

o 

0,0.0/ 

DATA  XZ1/ 

$  3.140,  3.261, 

2.919, 

2.115, 

1.531, 

0.960, 

0.615, 

0.400, 

0.230, 

$  2.809,  2.365, 

1.943, 

1.367, 

1.009, 

0.677, 

0.474, 

0.342, 

0.226, 

$  2.459,  1.665, 

1.254, 

0.776, 

0.510, 

0.286, 

0.164, 

0.096, 

0.047, 

$  2.070,  1.408, 

1.105, 

0.755, 

0.551, 

0.363, 

0.248, 

0.172, 

0.107, 

$  1.696,  1.142, 

0.910, 

0.624, 

0.447, 

0.280, 

0.181, 

0.118, 

0.067/ 

DATA  XZ2/ 

$  4.672,  4.300, 

3.695, 

2.532, 

1.755, 

1.053, 

0.656, 

0.419, 

0.237, 

$  4.465,  3.426, 

2.632, 

1.685, 

1.181, 

0.762, 

0.524, 

0.372, 

0.244, 

$  4.369,  2.726, 

1.909, 

1.093, 

0.694, 

0.380, 

0.216, 

0.124, 

0.061, 

$  4.297,  2.552, 

1.837, 

1.160, 

0.816, 

0.523, 

0.353, 

0.242, 

0.150, 

$  4.240,  2.468, 

1.805, 

1.147, 

0.798, 

0.490, 

0.314, 

0.203, 

0.115/ 

DATA  XZ3/ 

$10,090,  6.488, 

4.362, 

2.185, 

1.239, 

0.630, 

0.362, 

0.224, 

0.123, 

$  5.047,  2.653, 

1.604, 

0.812, 

0.525, 

0.328, 

0.223, 

0.157, 

0.102, 

$  3.095,  1.429, 

0.871, 

0.461, 

0.286, 

0.155, 

0.088, 

0.051, 

0.025, 

$  2.270,  1.127, 

0.771, 

0.478, 

0.336, 

0.215, 

0.146, 

0.100, 

0.062, 

$  1.983,  1.087, 

0.784, 

0.494, 

0.344, 

0.211, 

0.136, 

0.058, 

0.050/ 

WRITE (6, 100) 

100  FORMAT (*  SINGLE  EDGE  CRACK  PANEL  IN  TENSION  - * 

$  * PLANE  STRESS*,//) 

ITYPE=5 
I FLAG =1 
NUMAB=5 
NUMN=9 
DO  10  1=1,9 
X(I)-XX(I) 

X ( I ) =XY ( I ) 

10  CONTINUE 

DO  20  1=1,45 
Z1(I)=XZ1(I) 

Z2 (I) =XZ2 (I) 

Z3 (I) =XZ3 (I) 

20  CONTINUE 

RETURN 
END 
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APPENDIX  B 

A  SAMPLE  OUTPUT  LISTING  OF 
PROGRAM  EST 

(INPUT  IS  DEFINED  IN  APPENDIX  C) 


* 
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. .REWIND, TAPE6 
. .COPY ,TAPE6 
MATERIAL:  COPPER 

COMPACT  TENSION  -  PLANE  STRAIN 


o 

04 

u 


CM 

o 

I 

u 


u 

04 


u 

w 


p***-taoir»rn<N^ONaor*r*vfiir>^,,<rro<^cMrs 

aoao  p'l^r^r'p'vfivo'flvovo'fi'fivovovovflvo 

i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  1  i  1  1 

WWUWWWWCdWWUXCdXCdXMMX 

(NnmnHino'^uiNnhrnoHino^yfi 


^inifiinininininiPininiDinirinin^r 
ooooooooooooooooo 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

uuuutuuuuuuuuuuuuu 

r^minr>^OfNrr^r'(T>Hr)^vocDin 

HCJoo,»ror‘noMnH(^rrovo(NCDO 

vfiHHfNnn^vift\OvOhffl®ON^'-t 


o  o 

I  I 

u  u 

-h  r- 


yfiOr-^NHOMDr^vouMn^nn 
®®r*r»r»p^vovovovovovovovOvo 
i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 
CJ  MWWCJWWWCdMWWWXWCiJ 

Cu  Hno»/>0^(NonHoinooHH 

^,vflfl,vonoo(j\r*0'<,'flMi^oO 
fS<NOO(NrJOOHM(N(NHa>^,MOD 


^i^iniPininiruniPiniPin^Trv 

ooooooooooooooo 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
CJ  uuuuuuuuuuuuuuu 

U  ^•aNaoaor^r*ioinin*r^r<v'n>'-*0' 

fsHfsnn^invor^r'aoo'HHH 


< 

X 

04 

J 

< 


z 

04 


oooooooooooooooooooo 


z 

04 


oooooooooooooooo 


z 

x 


z 

x 


u 

J 


o 

o 

o 


oooooooooooooooooooo 
voinininminmininininminminin^,HT’*r 
ooooooooooooooooooo 
i  i  i  i  i  i  i  i  i  i  i  i  i  •  i  *  •  *  • 

uuuuuuuuuuuuuuuuuuu 

HCNoo^or-nounHr-vovccNCDOHH 

v0HHfSnr,X,^,invD^P'00®as3'HHH 


oooooooooooooooo 

cu  vaininirunirnmninininin-«r<,'r 

cn  ooooooooooooooo 

m  1  i  «  1  1  1  »  1  t  i  ■  1  1  •  1 

O  XXXXXXXXXXXXXXX 

j  ^•ONCDcorsr,*voiriinv^,mnHO' 

r*HfNjmn<«TirMafxhco^HHH 


X 

H 

Q 


1 ft 
cn 
o 


Q 

O 

X 

u 


ooooooooooooooooooo 
1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  •  »  1 

MUUUUCdUWUUMmJUUWUtdCC 

5r<n^r^^vooiri^^'v£)p'(7'OfNm5rvor^ 

f-4(NVinr'QDOHOHNn5TvO^CO(J'OH 

HCNrOTTlTUOCDONHHHHHr-IHHHMCM 


ininininininm,«r,*r^r'*rHrTT,*r^r 

ooooooooooooooo 
o  I  I  I  I  I  I  I  I  •  1  I  1  I  '  I 

o  UUUUUUUUUUUUUUU 

X  r^rcj^yD^oo^^Hinoooj^ 

Cj  nr‘Hvoo(N^H(Nroipvfir>*a'0 

HtN^invDOOCrHHHHHHHfN 


® 

o 

+ 

x 

o 

o 


u 

0 

0 

H 

X 

O 

eN 

0 

O 

cn 

04 

• 

r» 

o 

U 


x 


o 

+ 

x 

o 

o 


o 

o 


ao 

GO 


O 

O 

in  o 
cn  m 
•  in 
in 


X 

X 

Eh 

\ 

M 

<  It 

z 

3 

CN 

X 

X 

O 

X 

z 

O 

M 

< 

O 

4-3 

cu 

• 

0 

O 

cn 

TT 

ID 

Z 

ID 

x 

0  m 

0  • 

z 

• 

M 

O 

CN 

X 

Eh 

J 

O 

O 

X 

00 

X 

Eh 

ON 

X 

HT 

II 

X 

• 

II 

X 

Z  X 

X 

< 

Q  Q  Q  Q  Q  Q  Q  Q  o  H  H  H  iH  H  ^  N  (N  04  04  n 

noooooooooooooooooooo 
\oooooooooooooooooooo 
•^oooooooooooooooooooo 
^oooooooooooooooooooo 
^oooooooooooooooooooo 
^oooooooooooooooooooo 
X  ©  hhhhhhhhhhhhhhhhhhh 

_  HH(NNN{NN<NNNN(NN(N(N(SNN(N 

S  ooooooooooooooooooo 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
—  UUUUUUUUUUUWUWUWWUU 
ooinoinoinoinomoinoinoinoin 
inoo^^HrrooHinaDcsinc^oivoa\w^ 

nPK^HHfSN(imnn<rw»nini^\D^ 


ioinin^r^r^r^*',*»*,,5,,^,^^'rornrr,rnr^rr>ro 

ooooooooooooooooooo 
1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  <  1  1  > 

UUUUDQWWWUWUWWWUKlCdCCU 
vONOin^VDHOHvOirinH^COOOOnh 

^(NinHH(Sn<rinvOhO'^HHHf4(NM 


HH(NN(N(NfMfMfNCMMM(NCNWN(NMM 

ooooooooooooooooooo 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
uhuuuuuumwuwuwwuuuw 

ooinoinoinoinomoinoinoinoin 

mooir^HvcoHinooNinff'N'fiO'mys 

nP.HHH(N(N(Mnnnw^mirun^vo 


oooooooooooooooooooo 

oooooooooooooooooooo 

Qomomomoinomoinoinoinoinoin 

^  . . .  .  .  4  I . 

Oo  ii— tcNCNmm-^r^rininvoior-r^aooONON 

J 


oooooooooooooooooo© 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  •  • 

o  [dUUauWUWHWUWWWUWUBU 

CXi  u->4— i^cNCNcom^^rinininvoiDr^r-cocooi 

\  oomomomomomomomomom 

0*  l/-i^f_HfN<Njcnro*«r^,mmmmn*r^oocDO>>cT' 


Ol(T'CT'CT*(T*(7'<T'^C3r'OlCr'CJ'Cr'^O^G>G'ff'^5' 

XONONONONcnoNcnoiONOiONONONONONCNOioiONCi 


^oooooooooooooooooooo 
s.  mmmmmmmmmmmmmmmmmmmm 

<(Sf\KNM(N(NC4fMrJ(NN(M<NNP4NfN(NN(N 

aiaiaiaia>ocr»cT>oocr>ooc7>oooooo 


CN 

o 

CD 

ID 


£jQQOQOOOOr‘H'—^r‘H'-H'"Hn4(N(N 

OOOOOOOOOOOOOOOOO 

\0000000000000000 

OOOOOOOOOOOOOOOOO 

’—'OOOOOOOOOOOOOOOO 

HOOOOOOOOOOOOOOOO 

^OOOOOOOOOOOOOOOO 

o . 

HHC4fM(NNC4C4(N(N(N(N(N(N(N 

OOOOOOOOOOOOOOO 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
uuuuuuuuuuuuuuu 
ooooooooooooooo 

OOCNvDO^OOCNvDOrrCOCNiOO 
vr^HHNNNnrMf 


3 

CJ 

a : 

S 


II 

m 

X 


o 

O  cm 
m  m 

•  CD 


it 

X 
N 
<  II 

CN 

X 


o  m 
o  • 


o 

CN 


II 

fH 

Z  X 


®ininN,^J,^,Hr'*r^j,,3,c,r>roc''>mco 

ooooooooooooooo 
1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

uuuuuuuuuuuuuuu 

mcNor^^raoair-cninmminocN 

m^rr'cnrHO^H"*raiLnocNTr'sDai 

ConrsH(Mntrin^®HHHHH 


4— (»— ICNCNCMCNCNCNCNCNCNCNCNCNCN 
OOOOOOOOOOOOOOO 

XXXXXXXXXXXXXXX 

oaiooooooooooooo 

OO'NvDOVflONtfOVCDNtfO 


oooooooooooooooo 

oooooooooooooooo 

Qomomomomomomomom 

. . 

O  o  HHNNnnir^ininvovOM^ 


in 

M 

z 

o 

X 

cn 

>-* 

j 

o 

z 

CJ 


cn 

Eh 

J 

o 

cn 


x 

< 


o 

04 


ooooooooooooooo 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

uuuuuuuuuuuuuuu 

VOiaCOhVCvOlTlin'd’^COCMMH 

o>hm^ouTf- t  r*  cn  on  in  •— i  r**  cn  on 

LprHf-i(NCNmTr^rminmn‘r*cooD 


oionononononononononoicnonono>on 

XONONONONONOnOiONONONONONONONONON 

<mmmmmmmmmmmmmmmm 


^oooooooooooooooo 

\oooooooooooooooo 

<nnnnnpnnnnnnnnnnn 

OiONO'ONONONO'ONOiOiOiONO'OiONO' 

mmmmmmmmmmmmmmmm 


72 


CM  CM  ^ 
vO  v6  VO  VO 
till 

uuuu 
on  r*  m 
o>  r«*  o  o 
fN  ov  n  cr> 


ino'inco»-MO'®r'-®in^'r'»cocM.-4»-Hooo' 
aor'r‘M*'voovoo'Ovovo\ovovovoooin 
i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i 
U  UEdUCdUUUUCdUUUUUUUUUU 

x  r*n5rofNOHr'oor*0'0(Nrff>onr'n 

tN^H\OH(NvOr^OV(MN(NMr'OOVnH(N 

nfOHnnHNnnnrMH\ocN^n^trH 


vaD/ifNoovr*voir)«fl,nn(NHHooav(3v 
cor»r*r»r>ovoooovovoovoovooinm 
i  I  l  l  l  i  i  l  i  i  i  i  i  i  i  i  I  I  i 

(J  WWWWWWWWWWWWWWWWWWW 

X  ovr*iOiONfsr"r<TOM^Hnno^tx^i3v 

lOvoiTiMnoofScDOMnoHonncT'^roo' 

HiHlOHHlOHHHHHvO^HlTIHvOtNkO 


uuuu 
r-  in  co 
Nn^in 


inioimriirnfn/ufiinv-tfvvrv^vrr^^ 

0000000000000000000 
1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

U  WWWWWWWWWWWWWWWWWWW 

W  HOjm^invOf^OOOVHHHr-<HH(NCNfN(N 

oooooooooo^Hoar^^rmor^cDav 

HCNM^lOvOr^OOOVHHHHHHHHHH 


iflinininininmv^^^wvrww^ 

0000000000000000000 
1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 
CJ  uuuuuuuuuuuuuuuuuuu 

W  ®inco*HO'®^cMin®oco®o'CM^p»oco 

(Nin®Hn«ffvoH(N^invor-ONOHnv 
»cNcoin®r^®*H^Mf-irH,~i.--i^M.— icmcmcmcm 


Z 

X 


Z 

X 


OOOO 


OOOOOOOOOOOOOOOOOOOO 


OOOOOOOOOOOOOOOOOOOO 


z 

w 


z 

w 


OOOO 

till 

u  w  w  w 

MOnH 

(NPO^ifl 


OOOOOOOOOOOOOOOOOOOO 

Cu  inintOtnmLnLOtOin^r^r^r^r^r^rTr^r^r^r 

W  OOOOOOOOOOOOOOOOOOO 

tH  I  I  l  l  l  l  l  l  l  l  I  I  l  I  l  I  1  I  I 
Q  HUUUWUUCdUMUUUUUUUUCd 

J  HMnTTlDVDr'OavHHHHHHCNfNfMtN 

0000000000— icMcoTrinvop-®0' 
^OMCO^in®P^®ON'H^«-M-<^^M’H»-<'H— I 


OOOOOOOOOOOOOOOOOOOO 

cn  0000000000000000000 

»~i  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1  1 

Q  uuuuuuuuuuuuuuuuuuu 

J  QoinnH^v£>v(Nir®onvoO'fNvr^on 

(NinoHnvoavoHNvmvor^^OHnv 
H(Mr>in\Or*®HHHHHHHH(N(N(NfV| 


OOOO 

till 

uuuu 

O  ro  p*  H 
pvin^yfl 
CM  CM  CM  CM 


OOOOOOOOOOOOOOOOOOO 
Q  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

o  uuuuuuuuuuuuuuuuuuu 

X  ^QOCNvCOVinHOOM-Or^nOvCfSfTUDfvJ 

u  vOtNOMflNOHPnM'VOOOVHnM'Vfif^OVH 

HfOM>^®avHHHHHHrS<N(NNNCV»n 


irufiininiD^TrtfTrrrrfTrTrv^ 
OOOOOOOOOOOOOOO 
O  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

O  yuuuuuuuuuuuuuu 

X  p'-coov©co®®r^r'*r^®®®inm 

CJ  avONOvooHninr-ovHninp^ov 

Hninr^avHHHHHNCM(N(NM 


o  o 

I  I 

w  w 


o  o 

I  I 

w  w 


(Nponn 

OOOO 

OOOO 

OOOO 

OOOO 

OOOO 

OOOO 


OOOO 
OOOO 
o  m  o  in 


■  ® 
CM 


Sf 


CM  CM  CM  CM 
OOOO 
+  +  +  + 

2 

mm  w  w 

9'  O'  O'  O' 

H 

co  r*  h  m 

X 

®  ®  r— 

O 

•  •  •  • 

CM 

to 

O' 

® 

CO  CO  CO  CO 

CO 

OOOO 

• 

1  1  1  1 

CM 

W  W  W  W 
av  r*  r*  a> 

*h  ^  o 

CM  CM  CM  CO 
•  •  •  • 

•-3 

CM  CM  CM  CM 
OOOO 

n 

♦  +  +  + 

CO 

uuuu 

X 

ON  O  O'  O' 

co  h  in 

CO 

®  ®  r-  r- 

O 

WOOOOOOO— I— tCMCNOMCOCOCOTTTr 

•doooooooooooooooooooo 
\ooooooooooooooooo 00  o 

f-300000000000000000000 

—'OOOOOOOOOOOOOOOOOOOO 

H00000000000000000000 

0:00000000000000000000 

*-l»-lCNCMCMCMCMCMCMCMCMCMCMCMCNCMCMCMCM 
OOOOOOOOOOOOOOOOOOO 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
wwwwwwwwwwwwwwwwwww 
mav®<Mr^ro®-va,a\ino®<-Hr—  cm  i^no  v 

inom®<Nr-*-t®oino^,ONro®<Np-*H® 
^•a>«HHCMCMcncn^d,^*in»nin®®r-r^®® 


jninm^^^v^^nnnnnnnnnn 
OOOOOOOOOOOOOOOOOOO 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

uuuuuuuuuuuuuuuuuuu 

Hrri^r*r*corMoovcH^ravr>*r,*<Tinoo'av 

HM,a\r'r*<^M,o^Hnin»HM,®fMin^ 

H'VC7'H(Nnknr**CDHHHHfN(N(Nrr)rr)rrl 


— t«HCNCMCMCM<NCMCMCMCMCMCMCMCMCMCMCMCM 
OOOOOOOOOOOOOOOOOOO 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
uuuuuuuuuuuuuuuuuuu 

ino'v£>CMp,-cn®*ro'inovo— <p^<Nr-ro®^r 
mocomcMf'*— iiDoino^ONcoocMp^^H® 
-^^r-i^c^cMcofo^^rinininioior^t— ®® 


OOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOO 

Qoinomoinoinoinoinoinoinoinoin 

. . 

O  o  HHNCMfnnwmmvovor^r^ooovov 
J 


w 

6- 


r- 

co 

TT 

CN 


11 

co 

X 


o 

o  t- 

O' 


W©ooooo«-i^—»— <cMCMCMroco^r*Tinmso 

^OOOOOOOOOOOOOOOOOOOO 

\oooooooooooooooooooo 

^OOOOOOOOOOOOOOOOOOOO 
—'OOOOOOOOOOOOOOOOOOOO 
E-iOOOOOOOOOOOOOOOOOOOO 
XOOOOOOOOOOOOOOOOOOOO 
O* . 

WOHHHHHHHHHHHHHHHHHHH 

*-4CMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCMCM 
OOOOOOOOOOOOOOOOOOO 
W  +  +  +  +  +  +  +  +  +  +  +  +  +  +  ++  +  +  + 

uuuuuuuuuuuuuuuuuuu 

®^rinr— av<— tm^r®®o^rnmr-o>ofM^r 
^OinOin^VOHvOHr^(NMNMN®n® 
inHHcMfNnnTTTPininvcvor^r-ooovcr 


H 

X 

O 

w 


mm^v^^^vnnnnnnfnnnnn 
OOOOOOOOOOOOOOOOOOO 
I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
WWWWWWWWWWWWWWWWWWW 
^T^ravocrvr'^rav®'«T'«rr— m^Hmomin® 
vrr*<NnmHO'-i^vr(so^oNveHvoH 
HinHrMnin^c\HHHCNCMCNnnwin 


i-MCMCMCMCMCMCMCMCNCNCMCNCMCMCMCMCMCMCM 
OOOOOOOOOOOOOOOOOOO 
+  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  +  + 
WWWWWWWWWWWWWWWWWWW 
®**Tint—  O'— ico^rvoaoo— ic'Nini—O'OCM^r 
»HOinom«-t®*-t®—»r^cMr-cMr‘CM®ro® 
in— t— tcMCNcoco*ir««erinini£>®f^r-®®0'0' 


OOOOOOOOOOOOOOOOOOOO 

OOOOOOOOOOOOOOOOOOOO 

Qomoinoinoinomoinoinoinoinoin 

< . 

Oo  — i-H0M<Ncnco^<*rinin®®p-r^®®O'O' 

■J 


cn 

E- 


CO  CM  CM  CN 

\ 

►H 

^J*  CO  CO  CO  CO 

CO 

CO 

CO  CO  CO 

CO  CO 

CO 

CO  CM  CM 

CM 

CM 

CM 

\ 

M 

CO  CO  CO 

CO  CO 

CO 

CO 

CO 

CO 

CO 

CM  CM 

CM 

CM 

CM 

CM  CM 

CM 

OOOO 

<  N 

Z 

00000 

0 

0 

000 

0  0 

0 

OOO 

O 

O 

O 

<  n 

z 

•— s 

0 

OOO 

0  0 

0 

0 

0 

O 

0 

O  O 

O 

O 

O 

OOO 

till 

w  u  w  w 

CM 

X 

0 

u  u  u  u  u 

w 

1 

w 

WWW 

w  w 

w 

WWW 

w 

w 

w 

CM 

•mJ 

O 

W 

WWW 

w  w 

w 

W 

W 

W 

w 

w  w 

W 

w 

w 

w  w 

w 

HfH  P*  CO 

X 

X 

X 

P*  H  CM  CO  M* 

in 

vo  r- 

®  ON 

O' 

0  ®  CO 

0 

r» 

X 

X 

X 

CM 

O  ®  -H 

®  <— • 

r- 

CM 

p- 

CM 

CO 

CM  ^ 

ON 

00 

SO 

in  co 

CM 

in  o  o  — * 

w 

\ 

o  — t  ®  m 

CM 

ON 

iocoo 

r-  ^r 

— < 

O'  O  — t 

CM 

CM 

CO 

cn 

\ 

in 

r-*  m 

CM  »H 

O' 

® 

SO 

m 

CO 

O  *H 

»— 1 

CM 

m 

*«r  in 

SO 

ON  *"0  rH 

H-t 

X 

C*  H  (N  CM  n 

in  so  r— 

r-  ® 

O' 

O'  H  H 

r—i 

•— ( 

I— l 

*-• 

X 

® 

H  CM  CO 

in 

m 

SO 

P- 

® 

O' 

»-M  »-! 

H 

H 

f-4  »-H 

•—4 

•  •  •  • 

w 

CO 

0 

0 

0 

0 

r» 

z 

in 

z 

O'  O'  O'  O' 

CO 

w 

ON  O'  O'  O'  O'  O' 

O'  O' 

O'  ON  9V  O'  O' 

O' 

O'  O'  O'  O' 

O' 

O' 

® 

w 

ON  O' 

O'  O'  O' 

O'  O' 

ON 

O' 

O'  O' 

O' 

O'  O' 

O'  ON 

O'  O'  O' 

O' 

O'  ON  O'  ON 

O  CM 

W  O'  O'  O'  ON  O'  O' 

O' 

O'  O'  O'  O'  O'  O'  O' 

O'  O'  ^  O'  O' 

ON 

0  0 

WONO'O'ONONO'ONONO'O' 

O'  O' 

ON  ON 

O'  O' 

®  OS  Os  O' 

in  in  in  in 

O  • 

z 

<  VO  VO  ®  VO  V£>  ® 

VO  lO  vO  vO  VO  lO  VO 

®  VO  VO  vO  ®  ® 

v£> 

0  • 

z 

<  r-  r*  r—  t—  r- 

p-  p* 

p- 

f- 

r—  r—  r— 

r- 

r- 

r* 

f- 

•  •  •  • 

• 

M 

• 

w 

OOOO 
OOOO 
m  co  n  co 


O'  ov  ^  av 
O'  o'  O'  o' 
m  m  in  in 


o 

CM 


n 

H 
Z  X 


cn 

j 

x 

cn 

w 

x 

j 

j 

< 


^OOOOOOOOOOOOOOOOOOOO 
\mininininininininininininininmininin  in 
4!  CO  CO  Cl  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO  CO 


®VO®VO®VO®lOVO®VO®®®VO®VO®lD  SO 


o 

CM 


N 

fH 

Z  X 


cn 

E- 

j 

3 

£ 

W 

W 

< 


3:00000000000000000000 

\oooooooooooooooooooo 


ffv  iji  0V  0V  OV  9V  9v  O'  (Ji  OV  ON  9^  ^  9^ 

<C  9V  9V  9V  O'  OV  9v  O'  O'  O'  O'  9V  9v  9*  9V  9'  9'  9V  9V  9V  O' 

r-r^r^p-r'p'-r'r^r^r— r-r'C^r^p^r^p^r^r-r— 


73 


APPENDIX  C 

A  SAMPLE  INPUT  LISTING  OF 
THE  INTERACTIVE  COMPUTER 
PROGRAM  EST  FOR  A  MONOTONICALLY 
LOADED  CT  SPECIMEN 
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COMMAND-ATTACH , EST , ID=TUSIT , SN=AFML ,CY=1 
COMMAND-EST. 


ENTER  SPECIMEN  GEOMETRY 

(CCP ,CT  ,  DCP , STB , SET) == I CT 

PLANE  STRESS  (0)  OR  PLANE  STRAIN  (1)  ?1 

MATERIAL  == | COPPER 

E(NGLISH)  OR  M(ETRIC)  ?E 

IF  THE  FOLLOWING  VALUES  ARE  CORRECT, 

RESPOND  WITH  A  "Y"  .  IF  THEY  ARE  INCORRECT, 

ANSWER  "N",  THEN  ENTER  THE  CORRECT  VALUE 

MINIMUM  LOAD  (KIPS/N)  =  .0  (Y/N) ?Y 

LOAD  INCREMENT  (KIPS/N)  =  .25000  (Y/N) ?N 

0.5 

MAXIMUM  LOAD (KIPS/N)  =  5.3  (Y/N)?N 

10 

INITIAL  A/B  *  .25000  (Y/N)?Y 

RAMBERG -OSGOOD  EXPONENT  -  N  =  36.000  (Y/N) ?Y 

YIELD  STRENGTH  -  SIG0  (KSI/PA)  *  33000.  (Y/N)?Y 

PLASTIC  ZONE  CORRECTION  FACTOR  -  PZC=  1.0000  (Y/N)?Y 

ELASTIC  MODULUS  -  E  (KSI/PA)  =  .17000E+08  (Y/N) ?Y 

HALF  SPAN  (IN/M)  =  2.0000  (Y/N) ?Y 

WIDTH  -  W  (IN/M)  *  1.9950  (Y/N) ?Y 

THICKNESS  -  TH  (IN/M)  =  .49800  (Y/N) ?Y 

CRACK  LENGTH  INCREMENT  (IN/M)  =  .10000  (Y/N) ?Y 

CORR.  FACTOR  FOR  LIMIT  LOAD  -  LC  =  1.0000  (Y/N) ?Y 

MAXIMUM  DISPLACEMENT  (IN/M)  =  .10000E-01  (Y/N) ?N 

0.1 

STOP 

030400  MAXIMUM  EXECUTION  FL. 

.438  CP  SECONDS  EXECUTION  TIME. 
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